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SUMMARY 


The  topographic  accuracy  of  the  AN/APQ-102 (XA-2)  side-looking 
ra  ar  and  its  specific  applicability  to  1:50,000  and  1:250,000  scale 
topographic  mapping  were  tested  using  stereo  radar  techniques  with  real 
t  JW°  ^aSic  s^ereo  configurations,  opposite-side  and  same-side, 
es  e  un  er  various  conditions  relating  to  the  data  reduction 
echnique ,  control,  and  image  coordinate  weighting.  Also  tested  was 

e  geometric  fit  of  each  individual  radar  record  to  the  established 
ground  control. 

Although  original  plans  called  for  a  test  area  25  nautical 
miles  long  by  4  nautical  miles  wide,  shortcomings  in  the  data  actually 
acquired. led  to  adoption  of  a  much  smaller  test  area.  Moveover,  lack 
j  C?r  ain  co  lateral  data,  such  as  range  marks  and  time  marks,  required 
development  of  a  special  reduction  technique  based  upon  ground  control 
in  the  test  area. 


Photogrammetric  aerotriangulation  of  aerial  photography  flown 
simultaneously  with  the  radar  provided  a  means  of  defining  the  flight 
pa  Ground  points  in  the  test  area  were  also  aerotriangulated,  and 
used  as  control  and  as  the  standard  for  measuring  radar  accuracy. 

The  test  results  were  surprisingly  good,  with  the  opposite- 
side  stereo  configuration  providing  the  better  accuracy.  Average  root 
mean  square  (RMS)  values  of  7.7  meters,  12.1  meters  and  13.2  meters  in 
(across-track) ,  Y  (along-track)  and  Z  (elevation),  respectively,  were 
20  ^meters''  ^  oPPosite-side  case;  corresponding  RMS  values  of  9.5  meters 
™  !*rS  16-7.meters  were  obtained  for  the  same-side  case.  Results 
also  Hp  ?  the  individual  radar  records  to  the  aerotriangulated  positions 
demonstrated  excellent  geometric  fidelity  in  the  radar  imagery. 

.  Although  based  upon  limited  data,  the  test  results  indicate 
hat  the  AN/APQ-102 (XA-2)  radar,  using  stereo  techniques,  has  much  pro- 

or  all  weather,  day-night  mapping  at  scales  of  1:50,000  and  1:250,000. 
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FOREWORD 


This  research  was  undertaken  by  the  Autometric  Operation 
of  Raytheon  Company.  It  was  supported  by  the  Advanced  Research  Projects 
Agency  of  the  Department  of  Defense  under  ARPA  Order  No.  1229,  and  it 
was  monitored  by  the  U.S.  Army  Engineer  Topographic  Laboratories  under 
Contract  No.  DACA76-69-C-0002 ,  dated  12  March  1969. 

This  final  report  covers  the  period  from  March  1969  to  April  1970 
During  the  first  months  of  the  project,  the  Government  was  represented 
by  Mr.  Boyd  E.  Newman  of  the  Engineer  Topographic  Laboratories.  Later, 
the  duties  of  Government  Project  Engineer  were  assumed  by  Mr.  Frank 
Klimavicz,  also  of  the  Engineer  Topographic  Laboratories. 

The  project's  Principal  Investigator  was  Dr.  Gordon  Gracie  of 
Autometric/Raytheon.  Messrs.  Eldon  D.  Sewell  and  Harl  Hockeborn  acted 
as  in-house  project  consultants.  Other  Autometric/Raytheon  personnel 

Ta.  n.®  r®sP°nsible  Part  in  the  investigation  were  Messrs.  Ronald  K.  Brewer, 
John  W.  Bricker  and  Robert  A.  Johnson. 
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INTRODUCTION 


1.1 


Background 


has  Iona  k  evelopmeat  °f  an  all-weather,  day-night  mapping  capability 
ever  ,een  ^.desire^  Soal  of  the  military  mapping  community.  How¬ 

to  be  ’Jh™  ahlu  par^lcular  capability  may  yet  be,  it  is  most  likely 
leved  through  use  of  side-looking  airborne  radar  (SLAR).  One 

3r  SJstam’  th<f  APQ-97 (XE-1)  high-resolution  SLAR,  has  already  been 

Provinc  °  pY  applied  t0  1:250»000  scale  topographic  mapping  of  Darien 
Province,  Panama,  a  perpetually  cloudbound  region  of  the  earth. 

available  f ter,  sevaral  generations  of  SLAR  developr,  nt,  there  are  now 
available  a  number  of  coherent  side-looking  radar  systems  which  show 

ZaVqToITxA^  day“nl"ht  mappinR  ^-ations .  The 

/  Q  102 (XA  2)  radar  is  one  such  system,  and  its  application  to 

;  and  JarRa  Scale  topographic  mapping,  using  stereo  radar  tech¬ 
niques,  is  the  focal  point  of  this  investigation. 

.  I|lany  stuflles  have  already  been  undertaken  to  determine  the 
theoretical  feasibility  of  radar  mapping.  Radar  image  geometrv  and  its 
differences  from  that  of  the  central  projection,  error  propagation 
optimum  flight  line  geometry,  radar  parallax,  and  equipment  for  stereo 
radar  mensuration  are  some  of  the  topics  which  have  been  recentlv  in¬ 
vestigated  to  give  insight  into  the  radar  mapping  problem.  Little 
however,  has  been  done  in  the  way  of  practical  testing  of  real  radar 

^  dete™ine  actual  accuracy  and  develop  useable  techniques 

problem  whr  ,U  n^  “  thiS  particular  aaP-t  of  the  radar  mapping 
problem  which  is  addressed  in  this  investigation.  F  g 

.  Thl?  project  was  originally  designed  in  anticipation  of  pro- 

historvfatl°n"  FOll“'ln"  iS  3  brUf  rcsu"e  °f  Project's  Lrly 


Effective  date  of  contract  was  12  March  1969.  The  original 
project  plan  stipulated  SLAP.  Imagery  of  a  test  area  comprised  of  flat 
rolling  and  mountainous  terrain,  to  be  acquired  bv  the  Government  as  ’ 

photo35  P°SS1^e  after  c°ntract  initiation.  In  addition,  vertical  aerial 
Ph  tography,  flown  simultaneously  with  the  radar,  and  also  separately 
over  the  test  area,  was  to  be  provided.  This  photography  was  to  be 
taken  with  a  cartographic  camera,  and  some  means  was  to  be  provided  for 

AnriiC19r6r9  aT  tJe.alJulta"ao“a  photography  to  the  radar  record!  In 

as  Jnpit  for  th  ^  ^  th3t  ^  °f  thG  supporting  data  expected 

as  input  for  the  investigation  would  not  be  available.  Specifically 

there  would  be  no  time  marks  or  range  marks  on  the  radar  record  a^’no 
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correlating  data  between  the  cartographic  camera  exposures  and  radar 
record.  Autometric  was  also  notified  by  the  Government  that  the  test 
material  would  not  be  made  available  until  after  1  June.  The  project 
work  plan  was  revised  to  reflect  the  expected  changes. 

On  20  May,  the  Government  Project  Engineer  requested  Autometric 
to  stop  work  on  the  project  because  of  delays  in  data  acquisition  and 
uncertainty  as  to  the  nature  of  the  data  to  be  provided.  Autometric 
agreed  to  this  request,  and  no  further  work  was  done  until  18  September, 
when  the  radar  imagery  and  supporting  data  were  finally  received  from 
the  Government.  A  further  revision  of  the  project  work  plan  was  then 
made  to  accommodate  the  actual  data,  as  received,  and  the  contract 
expiration  date  was  extended  to  12  May  1970,  in  accordance  with  the  de¬ 
lay  experienced  in  receiving  the  data. 

The  incompleteness  of  the  acquired  data  made  it  quite  evident 
that  some  bootstrapping"  would  have  to  be  done  to  establish  the  missing 
information.  Although  it  was  expected  that  such  "bootstrapping"  would 
affect  the  results  (and  perhaps  their  validity)  to  some  extent,  it  was 
believed  that  useful  information  could  still  be  derived  from  the  inves- 
tigation.  The  actual  results  later  obtained  indeed  bore  this  out. 

1. 2  Objective 

The  objective  of  this  investigation  is  to  determine  the 
accuracy  with  which  topographic  map  information  can  be  extracted  from 
specific  stereo  imagery  obtained  from  the  AN/APO-102 (XA-2)  radar,  and  in 
particular  to  determine  the  applicability  of  the  AN/APO-102 (XA-2)  radar 
to  1:50,000  and  1:250,000  scale  topographic  mapping. 

1*3  Test  Material 


Radar  imagery  and  supporting  ..erial  photography  were  acquired 
by  the  4416th  Test  Squadron  of  the  USAF  Tactical  Air  Reconnaissance 
Center  (TARC) ,  located  at  Shaw  Air  Force  Base,  South  Carolina. 

The  or  ginal  flight  plan  specified  coverage  of  a  target  area 
25  nautical  miles  long  by  4  nautical  miles  wide,  in  the  locality  of 
Atlanta,  Ga.;  Chattanooga,  Tenn.;  or  Knoxville,  Tenn.;  the  final  area 
selection  being  dependent  on  weather  and  air  space  availability. 

Data  acquisition  was  to  be  conducted  in  two  parts  under  the 
following  specifications: 
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Tart  1:  High  Altitude 


Sensors : 

Altitude : 
Airspeed: 
SLAR  Mode: 
Offset : 
Passes : 


AN/APO-102 (XA-2)  Side-Looking  Radar 
KC-1B  Mapping  Camera  (stabilized) 

KS-72  3-inch  Vertical  Camera 

30,000  feet  Above  Ground  Level  (AOL) 

420  knots 

3  (optimum  side) 

6  nautical  miles  (NM)  from  near  edge  of  target  area. 

Two  passes  on  each  side  of  the  target  in  race  track 
pattern. 


Part  2:  Low  Altitude 
Sensors:  Same  as  Part  1 

5,000  feet.  Above  Ground  Level  (AGL) 
Same  as  Part  1 
Same  as  Part  1 


Altitude : 
Airspeed : 
SLAR  Mode: 
Offset • 


Passes : 


6  nautical  miles  (NM)  and  8  nautical  miles  (NM)  from 
near  edge  of  target  area. 

Two  passes  at  each  offset  distance  on  one  side  of 
the  target  area. 


was^ann^d’  ^  ph0t°8raPhic  run  5,000  feet  AGL  and  420  knots 

was  pianned  over  the  target  area  centerline,  using  onlv  the  KC-1B  and 
KS-72  cameras  (i.e.,  no  radar). 

The  flight  plan  was  amended  to  include  one  additional  mission 
a  pass  flam  at  30,000  feet  AGL  and  420  knots,  offset  18  NM  from  “he 
near  edge  o  the  target  area,  and  operating  in  SLAR  Mode  4. 

The  fli§ht  configuration,  as  planned,  is  shown  in  Figure  1 

degree)  ^  ^  *U  ^  P3rallel  (w±thin  ™ 

that  =Ld  cover  he  m  o“  l^ss  area,  and 

AM  „  ^The  area  actually  flown  for  the  test  was  the  one  located  near 

tint  '  ^e°fglac  33  Sh°Wn  in  FigUre  2*  Radar  and  aerial 

photography  for  Sorties  08-61,  08-75  and  08-212,  flown  on  8  August 

ves tigation?"  reSI’ectlvel>' •  for  L  in- 
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Generally  speaking,  the  flights  followed  the  flight  design 
given  in  Figures  1  and  2  fairly  well.  The  only  significant  departure 
from  the  plan  related  to  the  photography-only  flight  made  directly  over 
the  target  area.  This  pass  was  flown  at  high  altitude,  rather  than  at 
low  altitude  as  originally  planned,  probably  to  ensure  adequate  lateral 
coverage  of  the  target  area  with  the  KC-1B  camera. 

1*  4  Scope  of  the  Investigation 

This  .investigation  is  essentially  a  test  of  the  mapping 
accuracy  of  the  AN/APQ-102 (XA-2)  side-looking  radar,  using  real  data 
taken  under  real  flight  conditions  over  a  prescribed  target  area.  The 
test  not  only  includes  evaluation  of  planimetric  accuracy  but  also 
elevation  accuracy  as  well  under  stereo  radar  conditions.  The  two  basic 
stereo  radar  configurations  (opposite-side  and  same-side  cases)  are 
considered. 


An  essential  prerequisite  for  the  test  is  development  of  a 
suitable  procedure  tailored  to  handle  the  test  material  actually  ac¬ 
quired.  Development  of  the  test  procedure,  along  with  performance  of 
the  test  itself,  is  included  within  the  scope  of  the  investigation. 


The  project  is  divided  into  the  following  ten  tasks: 


Task 

1: 

Task 

2: 

Task 

3: 

Task 

4: 

Task 

5: 

Task 

6: 

Task 

7: 

Task 

8: 

Task 

9: 

Task 

10 

Inspection,  Evaluation  and  Selection  of  Test  Material. 
Design  of  Test  Procedure. 

Mathematical  Analysis  and  Programming. 

Radar  Mensuration. 

Aero triangulation  of  Photography  Flown  Simultaneously  With 
Radar. 

Aero triangulation  of  Photography  of  the  Test  Area. 

Radar  Data  Processing. 

Analysis  of  Results. 

Preparation  of  Interim  Technical  Report. 

Preparation  of  Final  Technical  Report. 


Tasks  1,  2  and  3  are  directed  toward  development  of  a  suitable 
test  procedure.  Task  1  assesses  the  quantity,  quality  and  completeness 
of  the  material  provided  for  the  test,  and  selects  a  specific  test  area; 
Task  2  designs  the  test  procedure  on  the  basis  of  information  found  in 

Task  1;  and  Task  3  implements  the  design  by  developing  the  necessarv 
software. 


Tasks  4  through  8  perform  the  test.  Tasks  4  through  7  are 
accomplished  in  accordance  with  the  test  procedure  developed  in  Tasks  1 
2  and  3;  and  Task  8  analyzes  the  test  results. 

Tasks  9  and  10  document  the  work.  The  Interim  Technical  Report 
was  prepared  in  February  1570,  following  some  initial  test  results.  The 
Final  Technical  Report,  this  report,  documents  the  entire  investigation 
including  the  work  described  in  the  Interim  Technical  Report. 
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INVESTIGATION 


^ 1  Inspection,  Evaluation  and  Selection  of  Test  Material  (Task  1) 

This  task  was  necessarily  the  first  to  be  accomplished,  since 
design  of  the  test  procedure  was  heavily  dependent  upon  the  quality, 
quantity  and  completeness  of  data  provided.  Upon  receipt  of  the  radar 
imagery,  vertical  photography  and  other  supporting  information,  the 
records  were  inspected  closely  and  organized  as  to  content  to  aid  in 
the  selection  of  specific  material  for  the  investigation . 

Initial  inspection  of  the  data  directed  attention  to  several 
problem  areas  which  might  at  least  hinder,  if  not  directly  affect,  the 
results  of  the  investigation.  Specifically,  for  the  KC-1B  vertical  frame 
photography  taken  simultaneously  with  radar  coverage  of  Sortie  08-61: 

a)  The  film  was  badly  fogged,  apparently  in  the  developing 
stage,  to  the  extent  that  much  resolution  was  lost, 
totally  in  some  areas . 

b)  No  camera  calibration  data  were  supplied  for  the  camera 
used  (KC-lB-Lens  Cone  No.  316).  As  a  result,  no  compen¬ 
sation  could  be  made  for  film  distortion,  and  a  fictitious 
set  of  master  fiducial  coordinates  would  have  to  be  de¬ 
duced  from  comparator  measurements  in  order  to  make 
possible  the  transformation  from  comparator  readings  to 
plate  coordinates. 

c)  Only  three  fiducial  marks  were  imaged  on  each  frame,  and, 
of  the  three,  one  was  a  partial  image  of  irregular  shape. 

d)  Cloud  cover  was  considerably  higher  than  the  normally 
accepted  limit. 

e)  Trees  were  in  complete  foliage,  creating  difficulties  in 
point  selection  and  identification. 

f)  Frame  spacing  was  irregular,  ranging  from  about  400  meters 
to  over  6000  meters  (Approximately  95%  endlap  to  60%  end- 
lap)  . 

g)  The  frame  clock  was  not  running. 

Radar  imagery  was  judged  to  vary  from  good  to  poor,  having  in¬ 
consistent  resolution  from  one  strip  to  another.  Most  resolution  problems 
were  associated  with  image  streaking,  image  fading,  or  overall  poor 
contrast.  The  absence  of  time  marks  and  range  marks  on  the  radar  imagery 
indicated  special  problems  which  would  have  to  be  dealt  with  in  the  design 
of  the  test  procedure. 
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A  specific  objective  of  Task  1  was  selection  of  a  test  area 
and  corresponding  radar  and  frame  camera  imagery  to  be  used  in  the  in¬ 
vestigation.  The  flight  plan  (Figures  1  and  2)  presented  several  possi¬ 
bilities  in  forming  radar  stereomodels  of  the  test  area.  All  possible 
configurations  were  examined,  and  one  radar  stereomodel  was  chosen  for 
each  of  the  two  general  cases  (opposite-side  and  same-side)  of  stereo 
radar  geometry. 

The  opposite-side  case  chosen  (See  Figure  3)  consists  of  two 
passes  offset  six  nautical  miles  to  either  side  of  the  target  area, 
both  at  an  altitude  of  approximately  32,000  feet  AGL  (Actual  flight  altitude 
was  higher  than  planned).  This  configuration  is  geometrically  the 
strongest  available,  giving  range  intersection  angles  of  about  105° 
throughout  the  target  area.  The  two  flights  comprising  the  same-side 
case  are  offset  six  and  eighteen  nautical  miles  to  one  side  of  the 
target  area  (Figure  3),  both  at  an  altitude  of  approximately  32,000 
feet  AGL  (again,  slightly  higher  than  planned).  This  combination  was 
thought  to  offer  the  best  visual  stereomodel  of  those  available.  The 
low  altitude  radar  imagery  was  dropped  from  consideration  because  of 
its  poor  quality. 

The  opposite-side  case  was  designated  Stereo  Radar  Model  I; 
the  same— side  case,  Stereo  Radar  Model  II.  The  specific  test  area 
selected  is  a  rectangular  strip,  approximately  5  miles  long  by  3  miles 
wide,  shown  in  Figure  4.  The  specific  radar  records  which  were  selected 
are : 

(1)  Sortie  08-61,  Pass  3,  Channel  A. 

(2)  Sortie  08-61,  Pass  4,  Channel  B. 

(3)  Sortie  08-61,  Pass  5,  Channel  A. 

For  sake  of  brevity,  these  records  will  he  referred  to  simply  as  Passes 
3,  4  or  5,  or  as  Strips  3,  4  or  5.  They  are  shown  in  Figure  3. 

Passes  3  and  4  make  up  Stereo  Radar  Model  I  (opposite— side 
case)  ;  Passes  3  and  5  make  up  Stereo  Radar  Model  II  (same— side  case)  . 

The  original  plan  anticipated  using  the  entire  target  area 
for  the  test.  However,  in  view  of  the  inconsistent  quality  of  photog¬ 
raphy  and  radar  imagery,  it  was  decided  to  confine  the  test  to  a  small 
area  where  the  radar  image  quality,  at  least,  was  as  good  as  possible. 

The  displacement  of  the  centerline  of  the  test  area  with 
respect  to  the  target  area,  very  apparent  in  Figure  4,  is  caused  by 
using  only  one  channel  of  each  pass.  Both  channels,  side  by  side,  cover 
a  swath  somewhat  wider  than  the  target  area. 

KC-1B  vertical  frame  photography,  taken  simultaneously  with 
the  selected  radar  passes,  was  indexed  to  USGS  Topographic  Maps.  In 
addition,  a  two-frame  model  of  vertical  photography  covering  the  test 
area  was  identified  and  selected  from  Sortie  08-212.  The  KS-72  frame 
photography  did  not  have  to  be  used. 
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FIGURE  3  STEREO  FLIGHT  CONFIGURATIONS  ACTUALLY  USED  IN  THE  TEST 


2.2 


The  original  test  approach  presumed  that  each  radar  record 
would  have  range  marks  in  the  across-f light  direction  and  time  marks  in 
ie  along-f light  direction  to  permit  direct  recovery  of  slant  range  and 
time  from  the  record  with  minimal  error.  Furthermore,  it  was  assumed 
that  precise  times  of  exposure  of  the  cartographic  camera,  operating 
simultaneously  with  the  radar,  would  be  suitably  indexed  on  the  corres¬ 
ponding  radar  record  to  permit  direct  correlation  between  photography 
and  radar  imagery.  Data  reduction  could  then  proceed  in  an  orderly 
manner;  that  is,  the  camera  air  station  positions  could  first  be  deter- 
mined,  then  the  radar  air  stations  could  be  interpolated  from  the  camera 
stations  ^th  aid  of  the  recorded  time  marks,  then  the  slant  ranges  could 
extracted  directly  from  the  radar  record,  and  finally,  the  ground 
positions  of  the  test  points  could  be  intersected. 

Since  the  actual  radar  imagery  which  was  received  had  neither 
range  marks,  nor  time  marks,  nor  any  index  marks  to  correlate  the  car  to¬ 
ff*?  camera  exposures  to  the  radar  imagery,  an  approach  different  from 
that  just  described  had  to  be  adopted.  Specifically,  it  was  necessary 
to  introduce  an  artificial  time  scale  and  make  use  of  ground  control 
ocated  in  the  test  area.  The  ground  control  points  have  known  ground 
positions,  and  are  also  imaged  on  the  radar  records.  This  provides  a 
tie  between  the  radar  records  and  the  ground,  in  much  the  same  wav  that 
lation  C°ntr01  tleS  aerlal  Photography  to  the  ground  in  an  aerotriangu- 

Inf ormation  derived  from  the  ground  control  can  be  used  to 
locate  the  radar  air  stations  along  the  flight  path  and  to  calibrate  the 
radar  range  measurements.  In  this  particular  approach,  the  camera  air 
stations  serve  only  to  define  the  flight  path,  independent  of  time.  In 
other  words,  without  any  correlation  between  camera  exposures  and  radar 
record,  radar  air  stations  cannot  be  time-interpolated  between  the  camera 
stations,  instead,  the  ground  control  is  used  to  locate  specific  radar 
air  stations  along  the  flight  path,  from  which  the  remaining  radar  air 
stations  are  interpolated. 

°nce  radar  air  station  positions  and  slant  ranges  have  been 
proach  lntersectlon  of  the  test  Points  proceeds  as  in  the  original  ap- 


The  specific  steps  taken  are  now  outlined: 

1.  Lay  off  a  reference  line  at  the  side  of  each  radar 
image  strip  (Figure  5).  This  line  can  be  defined 
by  two  finely  marked  holes  (reference  points)  in 
the  emulsion.  The  reference  line  serves  as  the 
t-axis  of  the  radar  image;  the  r-axis  is  defined  to 
be  perpendicular  to  the  t-axis  and  to  pass  through  the 
first  reference  point. 
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FLIGHT 


FIGURE  5.  MENSURATION  OF  RADAR  IMAGERY 


2.  Select  and  mark  a  number  of  points  which  are  identi¬ 
fiable  in  each  radar  image  strip  and  in  the  test  area 
photography.  Select  several  of  these  points  as  con¬ 
trol  points;  the  remaining  points  are  used  as  test 
points . 

3.  On  each  radar  image  strip,  measure  the  (x,y)  coordi¬ 
nates  of  all  points  (control  points,  test  points, 

and  the  points  defining  the  reference  line)  with  a  high- 
precision  comparator. 

4.  Transform  the  image  coordinates  of  all  points  to  the 
(r«0  plate  coordinate  system.  The  r-coord.inates  are 
uncalibrated  measures  of  slant  range;  the  t-coordinates 
serve  as  artificial  measures  of  time. 

5.  Aerotriangulate  the  photography  flown  simultaneously 
with  the  radar  to  obtain  the  camera  station  positions 
in  space. 

6.  Aerotriangulate  the  test  area  photography  to  establish 
the  ground  positions  of  the  control  points  and  test 
points  within  the  test  area. 

7.  Interpolate  the  triangulated  camera  station  positions 
to  obtain  the  flight  path  in  space  for  each  radar  pass. 

8.  Locate  the  zero-Doppler  line  for  each  control  point. 

Then  locate  the  position  of  the  radar  air  station  on 
the  flight  path  at  the  base  of  the  zero-Doppler  line. 
(Note:  for  a  specific  point  on  the  ground,  the  zero- 
Doppler  line  is  perpendicular  to  the  flight  path  and 
passes  through  the  point,  as  shown  in  Figure  6.  The 
radar  air  station  is  considered  to  be  located  where 
the  radar  actually  was  in  space  at  the  moment  it  was 
looking  at  the  point  along  the  zero-Doppler  line). 

9.  Determine  the  slant  range  for  each  control  point,  using 
the  known  ground  and  radar  air  station  coordinates. 

At  this  point,  in  addition  to  the  given  ground  coordi¬ 
nates,  the  following  data  are  known  for  each  control 
point : 

(a)  Radar  air  station  position. 

(b)  Slant  range. 

(c)  The  t-coordinate  ("time"),  and  the 
r-coordinate . 

10.  Calibrate  the  r-coordinates  of  the  control  points 
against  their  known  slant  ranges.  The  calibration 
curves,  one  for  each  flight,  can  be  represented 
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FIGURE  6.  ZERO-DOPPLER  CONDITION 


graphically,  or  analytically  by  a  polynomial  function. 
This,  in  effect,  substitutes  for  the  range  marks. 


11.  Construct  for  each  flight  an  ephemeris  of  radar  air 
station  positions  in  terms  of  the  artificial  time  scale. 

12.  Using  the  slant  range  calibration  curves,  determine  the 
slant  ranges  for  all  test  points  from  their  r-coordinates . 

13.  Interpolate  the  radar  air  station  position  for  each  test 
point,  using  the  t-coordinate  of  the  test  point  and  the 
ephemeris  constructed  in  Step  11. 

14.  Intersect  the  position  of  each  test  point,  using  the 
computed  radar  station  positions  and  slant  ranges  from 
each  of  the  two  flights  comprising  the  stereo  configura¬ 
tion.  Note  that  Steps  1  through  13  deal  with  each  flight 
separately.  Step  14  finally  merges  together  the  data 
from  both  flights. 

15.  Compare  the  intersected  positions  of  the  test  points  with 
their  positions  established  by  aerotriangulation  in  Step  6. 

Steps  1  through  4  are  the  essential  ingredients  of  Task  4 
(Radar  Mensuration).  Step  5  corresponds  with  Task  5  (Aerotriangulation  of 
Photography  Flown  Simultaneously  with  Radar),  and  Step  f>  corresponds  with 
Task  6  (Aerotriangulation  of  Photography  of  Test  Area).  Steps  7  through 
14  are  performed  in  Task  7  (Radar  Data  Processing),  and  Step  15  is  included 
in  Task  8  (Analysis  of  Results). 

2- 3  Mathematical  Analysis  and  Programming  (Task  3) 

The  salient  features  of  the  mathematical  analysis  and  program¬ 
ming  for  Task  7  (Radar  Data  Processing)  are  presented  in  this  section. 

More  detailed  mathematical  formulation  is  given  in  the  Appendix  of  this 
report. 


Section  2.3.1,  which  follows,  is  essentially  Step  7  of  the 
Test  Procedure;  Section  2.3.2  combines  Steps  8  and  9;  Section  2.3.3  is 

Step  10;  Section  2.3.4  is  Step  11;  and  Section  2.3.5  combines  Steps  12  13 
and  14.  " 

2-3.1  Interpolation  of  Flight  Path 


Camera  station  positions  along  the  flight  path  can  be  estab¬ 
lished  by  photogrammetric  aerotriangulation.  Other  points  on  the  flight 
path  can  then  be  located  by  interpolation  between  successive  camera  sta¬ 
tions.  The  camera  station  positions  and  interpolation  functions  thus  de¬ 
fine  the  flight  path  in  space. 
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Inspection  of  the  particular  test  material  to  be  used  in  this 
study  already  revealed  considerable  irregularity  in  camera  station  spacing 
(see  Section  2.1,  Item  (f)),  with  some  intervals  between  successive 
camera  stations  extending  more  than  6000  meters  (Figure  11,  appearing 
later  in  Section  2.5,  shows  the  irregular  camera  station  spacing).  It 
was  felt  that  the  excessive  separation  between  camera  stations  did  not 
justify  use  of  anything  other  than  linear  interpolation  between  successive 
stations.  Thus,  the  flight  path  was  defined  by  a  set  of  line  segments 
connecting  successive  camera  stations. 

2,3,2  Determination  of  Radar  Air  Station  Position  and  Slant  Range 

for  Each  Control  Point 


It  is  assumed  that  the  cartographic  camera  and  radar  antenna 
are  located  at  essentially  the  same  position  in  the  aircraft;  that  is, 
there  is  negligible  physical  displacement  between  the  camera  and  antenna 
center.  Hence,  the  flight  path  defined  by  the  camera  stations  is  also 
assumed  to  be  the  flight  path  of  the  radar. 

Once  the  radar  flight  path  has  been  defined,  the  radar  air 
station  for  each  control  point  can  be  located.  This  is  done  by  construc¬ 
ting  the  control  point's  zero-Doppler  line;  that  is,  a  line  perpendicular 
to  the  flight  path  and  passing  through  the  control  point  (Figure  6).  The 
base  of  the  perpendicular  locates  the  position  of  the  radar,  the  radar 
air  station,  at  the  time  the  control  point  was  imaged  on  the  radar  record. 
(Actually,  the  radar  imaging  process  for  a  coherent  system  is  more  in¬ 
volved  than  what  this  simple  geometrical  model  indicates  in  that  the 
radar  is  really  collecting  an  extended  phase  history  of  the  control  point 
as  the  aircraft  travels  over  a  sizeable  portion  of  the  flight  path.  How¬ 
ever,  subsequent  construction  of  the  radar  image  in  a  correlator-processor 
is  done  in  such  a  way  that  the  control  point  can  be  taken  as  viewed  simplv 
along  the  zero-Doppler  line.) 


_  The  radar  air  station  position  (X  Y  Z)  is  obtained  bv 

solving  the  following  three  equations:  K  K 


“(XR  -  XI>  +  «Yr  -  V  +  Y<Zr 

R 


V 


=  0 


XR  =  Xo  +  f  <Y»  -  Y-> 
-R  =  7  A 


z«  =  zo  +  F  <yr  -  V 


(1) 

(2) 

(3) 


in  which  a,  fc>  and  y  are  the  direction  numbers  of 
in  which  the  radar  air  station  lies;  X  ,  Y  and 
of  the  camera  station  which  initiates  the  flight 
Xj.  Yj  and  Zj  are  the  ground  coordinates  of  the 


the  flight  path  segment 
Z0  are  the  coordinates 
path  segment;  and 
control  point. 
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If  the  zero-Doppler  line  is  t_aken  normal  to  the  ground  track 
t f  the  flight  path,  rather  than  to  the  flight  path  itself,  the  third  term 
of  Equation  (1)  is  eliminated  first  before  solving  for  the  radar  air 
station  position.  For  a  nearly  level  flight  path,  there  should  be  little 
difference  between  using  the  flight  path  or  ground  track,  because  for  this 
situation  y  is  already  very  small. 

Once  the  position  of  the  radar  air  station  is  known,  the  slant 
range  (SR^)  to  the  control  point  can  be  computed: 


SRI  =  ^R 


V2  +  (yr  -  V2  +  <ZR  -  V2J 


(4) 


Slant  ranges  computed  for  several  control  points  will  be  used 
to  determine  the  slant  range  calibration  function. 


2.3.3 


Slant  Range  Calibration 


It  is  assumed  that  a  preliminary  transformation  (Step  4  of  the 
test  procedure  outlined  in  Section  2.2)  has  first  been  made  to  put  all 
image  coordinates  measured  in  the  comparator  (x,y)  coordinate  system  into 
the  plate  (r,t)  coordinate  sv.tem  defined  by  the  reference  points  (Figure 
5).  If  the  reference  points  have  been  cai'efully  chosen,  the  r-coordinate 
will  be  essentially  an  uncalibrated  measure  of  slant  range.  However,  to 
account  for  any  error  in  the  direction  of  the  selected  reference  line,  a 
new  coordinate  system  (R, T)  is  introduced,  Figure  7,  in  which  the  T-axis 
makes  a  small  angle  0  with  the  t-axis.  All  points  on  the  T-axis  are  pre¬ 
sumed  to  have  the  same  slant  range,  and  the  R-coordinates  can  then  be 
directly  related  to  slant  range.  Finding  0 ,  of  course,  is  an  added  re¬ 
quirement  . 


re 


The  plate  coordinates,  Rj  and  T  of  the  image  of  point  I  are 
lated  to  the  plate  coordinates  rr  and  t  by  the  following  equations: 

Rj  =  rl  cos  3  +  t  sin  0  (5) 


TI  -  tI  cos  ®  “  ri  sin  0 


(6) 


Now  the  slant  range,  SR  ,  for  point  I  is  assumed  to  be  a 
second-order  polynomial  function  o t  the  Rj-coordinate ;  that  is, 


SRT  =  A 


SR 


+  B 


SR  7. 


R,  +  C 


SR  RI 


(7) 


Tnls  is  an  empirical  relationship,  but,  nevertheless,  a  realistic  one. 

The  constant  term,  Ag^,  is  the  slant  range  for  points  whose  images  lie  right 
on  the  T-axis;  the  coefficient  Bg^  is  a  scale  factor,  relating  distances 
measured  along  the  R-axis  to  corresponding  slant  ranges;  and  the  coefficient 
CgR  accounts  for  any  linear  variation,  hopefully  small,  in  the  scale  across 
the  radar  record. 


/*  -AXIS 


T-  AXIS 

t- AXIS 


FIGUHE  7.  PLATE  COORDINATE  SYSTEMS  FOR  RADAR  IMAGERY 


i  .  -i  ^ant  ran£e  calibration  can  be  accomplished  by  fitting  the 
oolyncnial  (Equation  (7))  to  data  obtained  for  control  points  imaged  in 

h.  .a  ar  The  flt  ls  best  made  usang  the  method  of  least  squares, 

,  in  addition  to  the  polynomial  coefficients  A™,  B  and  C  ,  also 

editions"  d  Unkn°Wn,rotation  an«le.  6  •  The  for^of  SR  the  ^condition 
the  ri  Kt  1?Jthe  least  sc5uares  solution  is  obtained  by  substituting 

the  right-hand  side  of  Equation  (5)  for  Rj  in  Equation  (7);  thus, 

SRI  =  ASR  +  BSR^rl  cos  0  +  ti  sin  0)  +  cSR(rT  cos  6  +  tj  sin  0)2  (8)  . 

have  heen  ^  Parameters  (ASR,  BSR,  CSR  and  6)  of  the  calibration 

ItlJ  n ?  Y  5  Particuiar  ™dar  ^cord,  slant  ranges  for  all 

other  points  imaged  and  measured  (known  rT  and  t  )  in  the  record  can  be 

°  ained  using  Equation  8.  In  particular,  the  slant  ranges  for  the  test 
be  ^PUted.  However,  this  computation  is  not  actually  made 
this  time.  Instead,  it  is  incorporated  into  the  least  squares  inter¬ 
section  computation  which  follows  later. 

p  If  53dar  exbibitfi  stability  in  the  display  of  slant  range 

to  the  ^  CSR  remaln  alm°St  conatant  from  one  record 

the  next),  information  From  ‘  two  or  more  records  can  be  combined  to 

determine  common  values  for  BSR  and  C  .  The  other  parameters  A  and  8 

must  still  be  evaluated  for  each  indi^dual  record,  however,  sinJRthey  ’ 

depend  upon  selection  of  the  reference  points,  which  is  somewhat  arbitrary. 


2.3.4 


"Time"  Correlation 


A  necessary  step  in  the  data  reduction  is  correlation  of  the 
radar  air  station  positions  with  the  radar  record.  More  specifically 

station"55  "-J-  ^  °btained  which  can  be  used  to  determine  the  radar  air 
i  l  5  measured  image  coordinates  on  the  radar  record. 

Radar  air  station  positions  have  already  been  determined  for  the  control 

the"rada  ^  P°intB  have  been  identified  and  measured  on 

and  the  r  ^  '  <’  *  V°sitive  relationship  between  the  radar  record 

he  radar  air  stations  can  be  established  through  the  ground  control. 

Tbe  T_axi-S  defined  in  Section  2.3.3  can  be  regarded  as  an 
improved  time  axis,  m  the  sense  that  it  is  an  improvement  ovei  the 
t-axis  defined  by  the  arbitrarily  selected  reference  points.  The 
coordinates  are  still  artificial  measures  of  time,  though,  since  no  real 
time  was  actually  recorded  during  operation  of  the  radar.  However  it  is 
really  immaterial  whether  the  T-coordinates  are  real  or  artificial 
measures  of  time,  since  all  that  is  needed  is  an  independent  variable  for 
in  eqpoiauo,,  of  the  radar  air  stations.  Convenient!?,  the  pl«. loordl 
nate,  Tj,  la  selected  as  the  Independent  variable,  and  interpolation  of 

lTyl^llTt'°n  C°°rdi"ateS  -i".  the  following  -tS? 


2 


+ 


(9) 


Ax  +  Bx  Ti  +  cx  Ti 

\  =  Ay  +  BY  TI  +  CY  TI2  +  •••  (10) 


Az  +  Bz 


TI  +  CZ 


(11) 


in  which  Y^  ,  Z^  are  the  coordinates  of  the  radar  air  station  for 

the  I th  ground  point;  and  Ax,  Bx>  Cx,  Ay,  B  Cy,  Az,  Bz,  Cz,  etc,,  are 
polynomial  coefficients  which  have  to  be  first  evaluated. 


Lvaluation  of  the  polynomial  coefficients  can  be  accomplished 
by  fitting  the  time  polynomials  to  known  data  (Xp  ,  Yp  ,  Zp  and  T  )  for 
the  control  points.  A  separate  set  of  I  KI  RI  1 

coefficients  is  found  for  each  flight  segment,  using  the  method  of  least 
squares.  Because  each  flight  segment  is  a  straight  line,  the  A-coef f icients 
are  necessarily  related  in  accordance  with  Equations  (2)  and  (3),  as  are 
the  B-coefficients,  C-coeff icients ,  and  so  on.  Furthermore,  if  the  ratio 
of  aircraft  ground  speed  to  radar  data  film  speed  is  constant,  second  and 
higher  order  terms  in  the  "time"  polynomials  are  zero,  and  the  interpola¬ 
tion  of  radar  air  stations  along  the  flight  path  is  linear  with  respect 
to  Tj , 


Once,  the  polynomial  coefficients  have  been  found  for  each 
flight  line  segment,  the  radar  air  station  position  for  any  test  point 
can  be  computed  using  the  time  coordinate,  TT,  of  the  point's  image  on  the 
radar  record.  However,  as  with  the  point's  slant  range,  this  computation 
is  not  actually  made  at  this  time,  but  is  incorporated  into  the  least 
squares  intersection  computation  which  follows. 

2.3.5  Intersection  Geometry 

If,  for  any  given  test  point,  the  direction  of  flight,  slant 
range  and  radar  air  station  position  are  known,  or  readily  determinable 
for  each  of  two  radar  flights  making  up  a  stereo  radar  model,  the  position 
of  the  test  point  can  be  intersected  in  space.  Computation  of  all  data 
necessary  for  the  space  intersection  has  already  been  dealt  with  in 

Sections  2.3.1  through  2.3.4,  and  so  these  data  are  now  presumed  to  be 
available. 


For  any  test  point,  I,  imaged  in  the  two  radar  records  making 
up  the  stereo  radar  model,  four  equations  can  be  written  expressing  slant 
range  and  zero-Doppler  conditions  for  the  point's  intersection.  Specif¬ 
ically,  for  the  first  radar  record: 
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(12) 


[(xi  ■  v2  +  (Yi '  yrx)2  +  (zi  ■  \)2]h  =  srj1 

ai(XI  ■  +  ei(Yi  -  YRi>  +  Yx  (Z!  -  zRj)  =  0  (13)  , 

and  for  the  second  radar  record: 

t(XI  "  V2  +  (YI  "  V2  +  (ZI  -  \)2]h  =  SRj2  (14) 

a2(XI  "  +  32(YI  "  Yr2>  +  Y2(Zj  -  ZR^)  =  0  (15)  . 

In  these  equations,  X  ,  Y  and  Z  are  -he  coordinates  of  the  first 

Ki  R  l  Ri 

radar  air  station  (pertaining  to  the  first  radar  record),  and  Xp  ,  Y 
3nd  ZR2  are  the  c°°rdirates  of  the  second  radar  air  station*  ?V2  ’  Rz 

»  3i _  and  Yi  are  the  direction  numbers  of  the  first  flight  path 
(pertaining  to  the  first  radar  record),  and  a2  ,  R  and  y,  are  the 
direction  numbers  of  the  second  flight  path} 

SRj^  is  the  slant  range  to  I  from  the  first  radar  air  station,  and  SR-r 

i2 

is  the  slant  range  to  I  from  the  second  radar  air  station;  and  XT,  Y 
and  Zj  are  the  unknown  ground  coordinates  of  I .  11 

Equations  (12)  and  (14)  express  the  slant  ranges  between 
radar  air  station  R!  and  test  point  I,  and  radar  air  station  R2  and  I, 
respectively,  in  terms  of  the  object  space  coordinates  of  the  points. 
Equations  (13)  and  (15)  express  the  orthogonality  of  the  flight  paths  and 
the  corresponding  zero-Doppler  lines  in  terms  of  direction  numbers. 

Since  four  equations  are  used  to  solve  for  the  three  unknown 
coordinates  (Xj ,  Y_,  Zj) ,  a  least  squares  solution  is  employed.  In  this 
solution,  the  plate  coordinates,  Rj  and  Tj,  implicit  in  the  slant  range 
and  in  the  radar  air  station  coordinates,  are  taken  as  the  observations. 
Equations  (7),  (9),  (10)  and  (11)  are  thus  included  with  Equations  (12) 
through  (15)  to  form  the  condition  equations  of  the  least  squares  solu¬ 
tion.  By  setting  up  the  condition  equations  in  this  manner,  prior  compu¬ 
tation  of  slant  range  and  radar  air  station  coordinates  is  not  required. 

It  might  be  noted  at  this  point  that  two  positions  of  inter- 

b!  f°r  eaCh  t6St  P°int*  since  eacd  pair  of  equations 

((  2)/(13)  and  (14) / (15) )  def ines  a  circle  in  space,  and  the  two  circles 
thus  defined  by  all  four  equations  can  intersect  at  two  points.  In 
practice,  this  ambiguity  presents  no  problem  so  long  as  reasonably  good 

estimates  of  Xj,  Yj.  and  Z  are  given  in  the  first  place  to  initiate  the 
least  squares  solution. 
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2.3.6 


Error  Propagation 


f  A  useful  feature  of  the  least  squares  solution  is  provision 

model  for  Propagating  errors  in  the  observed  quantities  into  the 
.  nal  results.  The  error  propagation  model  is  derived  directly  from  the 
rse  o  the  coefficient  matrix  of  the  normal  equations  constructed 
urmg  the  course  of  the  least  squares  solution.  In  the  case  at  hand, 
measures  of  precision  for  the  computed  ground  coordinates,  XT ,  i'  Z  , 
can  be  obtained  from  the  least  squares  intersection,  based  upon  precisions 
ot  the  image  coordinates,  RT  and  Tt. 


The  measures  of  precision  for 
dependent  upon  the  precisions  assigned  R-r 
precisions  should  be  as  realistic  as  poss 
results. 


^1  *  Yj  and  Z  are  necessarily 
and  Tj .  Hence,  the  assigned 
ibis  to  make  any  sense  out  of  the 


Th(7  least  squares  intersection  solution  assumes  that  the  slant 
range  calibration  constants  and  other  determined  coefficients  included  in 
the  condition  equations  are  errorless,  which  is  not  entirely  correct. 

*Mpr^eleSS’/Vuit:<f  °f  tHiS  shortconin*.  error  propagation  is  still 
bL u  ,glve  insi8ht  into  the  precision  with  which  test  point  positions 
can  be  determined  using  this  particular  stereo  radar  technique! 


2.3.7 


Computer  Programs 


data 


Four  computer  programs  were  prepared  to  accomplish 
reduction.  Brief  descriptions  of  these  programs  are  given 


the  radar 


Program  1: 


This  computer  program  defines  the  flight  path,  and  deter¬ 
mines  the  radar  air  station  positions  for  the  known  control 
points  using  the  zero-Doppler  condition.  It  also  computes 
the  slant  ranges  for  the  control  points.  The  mathematical 
analysis  for  the  program  is  described  in  Sections  2.3.1 
and  2.3.2.  The  program  has  the  option  of  defining  the 
zero-Doppler  line  normal  to  the  flight  path  or  normal  to 
the  ground  track. 


Program  2; 


Program  3: 


tHe  alant  _  range  calibration  program.  In  accordancf 
w  tn  the  mathematical  analysis  given  in  Section  2.3.3 
this  program  uses  the  method  of  least  squares  to  find' the 
slant  range  calibration  parameters.  The  program  allows  for 
processing  all  radar  records  separately  or  together.  Com¬ 
putation  of  slant  range  residuals  is  included. 


This  program  performs  the  "time" 
Section  2.3.4.  Specifically,  it 
estimates  of  the  time  polynomial 
point  data. 


correlation  described  in 
computes  least  squares 
coefficients  from  control 


Program  4: 


This  is  the  intersection  program, 
analysis  outlined  in  Section  2.3.5, 


In  accordance  with  the 
this  program  computes, 
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by  intersection,  the  positions  of  the  test  points  on  the 
ground.  The  program  has  the  option  of  defining  the  zero- 
Doppler  line  normal  to  the  flight  path  or  normal  to  the 
ground  track.  A  least  squares  solution,  including  error 
propagation  (Section  2.3.6),  is  used. 

All  programs  were  written  in  FORTRAN  IV  for  the  IBM  360/44  computer. 

2*4  Radar  Mensuration  (Task  4) 

The  sections  of  radar  Strips  3,  4  and  5  of  Sortie  08-61,  mak¬ 
ing  up  the  two  models  under  consideration  (Model  I,  opposite-side  case. 
Strips  3  and  4;  Model  II,  same-side  case,  Strips  3  and  5),  were  contact 
printed  onto  glass  diapositive  plates.  Under  optical  magnification,  a 
number  of  points  (to  be  used  as  control  and  test  points)  were  carefully 
selected  and  marked  on  each  plate  for  measurement  on  the  Mann  comparator. 
All  points  selected  were  also  identifiable  on  the  aerial  photography  of  the 
test  area.  A  total  of  60  points  were  selected  on  Strip  3,  55  on  Strip  4, 
and  30  on  Strip  5.  Fewer  points  were  selected  on  Strip  5  because  there 
was  less  detail  appearing  on  that  strip.  Of  the  points  selected,  49  were 
common  to  Strips  3  and  4,  and  so  could  be  used  for  Model  I  (Figure  8). 

All  30  points  selected  on  Strip  5  also  appeared  on  Strip  3,  and  so  could  be 
used  for  Model  II  (Figure  9).  Twenty-five  points  were  common  to  both  radar 
models. 


Each  image  selected  for  measurement  was  marked  on  the  Wild  PUG 
point  transfer  device  by  drilling  a  50-micrometer  (ym)  hole  in  the  plate 
emulsion.  Radar  Strips  3  and  5  (Model  II,  same-side  case)  were  viewed 
stereoscopically ,  and  22  of  the  30  points  selected  on  Strip  5  were  drilled 
on  that  strip  in  this  mode.  Extensive  efforts  were  made  to  view  and  mark 
points  stereoscopically  in  the  opposite-side  case  (Model  I),  with  unsatis¬ 
factory  results.  Consequently,  the  points  on  Strips  3  and  4  were  marked 
monoscopically . 

A  reference  line,  was  established  on  each  radar  record  by  drill¬ 
ing  two  50-micrometer  (ym)  holes  in  the  plate  emulsion  at  each  end  of 
the  test  area,  along  the  near  edge  of  the  imagery.  This  reference  line 
constructed  as  nearly  parallel  to  the  edge  of  the  imagery  as  visual  means 
allowed,  forms  the  basis  of  the  radar  plate  coordinate  system  into  which 
the  comparator  measurements  are  transformed  (Figure  5).  The  along-flight 
direction  (the  reference  line  itself)  defines  the  artificial  time  scale  by 
which  image  points  can  be  correlated  with  radar  air  station  positions.  The 
axis  normal  to  the  time  axis  is  taken  as  the  range  scale. 

Coordinates  (x,y)  of  all  radar  image  points  and  reference  points 
were  carefully  measured  on  the  glass  diapositive  plates  with  the  Mann  com¬ 
parator  The  measurements  were  then  transformed  to  the  plate  coordinate 
system  (r,t)  defined  by  the  reference  points. 
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Northing  in  Maters 


Y 


FIGURE  8.  SELECTED  GROUND  POINTS-MODEL  I 
(OPPOSITE-SIDE  CASE) 
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731000 


Northing  in  Meters 


Y 


FIGURE  9.  SELECTED  GROUND  POINTS-MODEL  II 
(SAME-SIDE  CASE) 
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“ *  ^  Aerotriangulation  of  Photography  Flown  Simultaneously  With 

the  Radar  (Task  5) 

Vertical  photography,  flown  simultaneously  with  the  radar 
coverage  of  the  test  area,  was  photogrammetrically  aerotriangulated  to 
establish  the  flight  p/.th  of  the  aircraft  (and  of  the  radar)  at  the  time 
of  flying  the  radar.  Camera  station  coordinates  were  determined  for  each 
frame  of  photography  in  the  three  strips  selected  for  the  study  (Strips 
3,  4  and  5  of  Sortie  08-61),  giving  for  each  strip  a  set  of  discrete  points 
along  its  flight  path. 

Ten  frames  in  Strip  3,  9  frames  in  Strip  4,  and  14  frames  in 
Strip  5  were  selected  for  aerotriangulat ion.  The  frames  were  contact 
printed  onto  film  positives,  and  several  points  (9—21  per  frame)  were  se¬ 
lected  for  mensuration  and  marked.  Mensuration  was  then  accomplished  on 
a  Mann  precision  comparator. 

Preprocessing  of  image  coordinates  for  the  aerotriangulation 
had  to  be  tailored  to  accommodate  the  available  camera  data,  incomplete  as 
it  was.  No  compensation  could  be  made  for  film  distortion,  and  a  set  of 
master  fiducial  mark  coordinates  had  to  be  artifically  deduced  from  compara¬ 
tor  measurements  of  the  three  available  fiducial  marks  in  order  to  define 
the  position  of  the  principal  point  (Figure  10). 

A  lens  distortion  table  for  the  camera  (KC-lB,  Lens  Cone 
No.  316)  was  extracted  from  U.S.  Army  Geodesy,  Intelligence  and  Mapping 
Research  and  Development  Agency  Technical  Note  66-2,  "The  Calibration  of 
Military  Cartographic  Cameras",  and  applied  in  the  preprocessing  of  plate 
measurements.  The  corresponding  calibrated  focal  length  was  also  obtained 
from  the  same  source.  In  addition,  plate  coordinates  were  corrected  for 
effects  of  atmospheric  refraction. 

Ground  control  used  for  the  aerotriangulation  of  Strips  3  and 
4  was  scaled  from  1:24,000  USGS  quadrangle  maps  (UTM  grid),  and  that  for 
3  was  scaled  from  USGS  maps  at  1:62,500  scale.  Standard  deviations 
of  12  meters  in  X  and  Y  and  2  meters  in  elevation  were  assigned  to  the  con¬ 
trol  for  weighting  purposes. 

The  strips  were  aerotriangulated  using  a  well-established 
analytical  least  squares  solution.  Data  reduction,  including  error  propa¬ 
gation,  was  performed  using  Autometric's  in-house  IBM  360/44  computer. 

The  a  posteriori  variance  of  unit  weight  computed  from  the  least  squares 
solution  indicated  a  standard  deviation  of  0.040  mm  in  the  measurement  of 
each  image  coordinate.  This  seems  to  be  a  reasonable  value,  considering 
the  quality  of  the  photography  and  the  incomplete  camera  calibration  data. 

The  aerotriangulation  is  briefly  summarized  in  Table  I.  The 
triangulated  camera  station  positions  are  plotted  in  Figure  11. 
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TABLE  I.  SUMMARY  OF  AEROTRIANGULATEb  STRIPS 
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FIGURE  11.  CAMERA  STATION  POSITIONS 


2.6 


Aerotriangulation  of  Photography  of  Test  Area  (Task  6) 


A  single  model  was  selected  from  Sortie  OS-212  photography 
(Frame  2334,  seventh  exposure;  and  Frame  2336)  for  photogrammetric  aero- 
triangulation  of  test  and  control  points  in  the  test  area.  The  two  frames 
were  contact  printed  onto  glass  diapositive  plates.  The  images  of  all 
points  selected  on  the  radar  strips  in  Task  4  were  then  drilled  stereo- 
scopically  on  the  plates  with  the  Wild  PUG  and  measured  on  the  Mann  compara¬ 
tor.  Additional  points  lying  outside  the  test  area,  tn,t  within  the  con¬ 
fines  of  the  photogrammetric  model,  were  selected,  drilled,  and  measured 
in  order  to  strengthen  the  model's  geometric  solution.  A  total  of  74 
image  points,  plus  the  3  fiducial  marks,  were  measured  on  each  plate. 

Preprocessing  of  image  coordinates  was  accomplished  in  accor¬ 
dance  with  the  approach  outlined  in  Section  2.5.  Horizontal  and  vercical 
control  was  scaled  from  two  USGS  quadrangle  maps  (1:24,000  scale).  In 
addition  to  points  chosen  specifically  to  control  the  aerotri angulation, 
all  points  selected  in  Task  4  which  were  identifiable  on  the  maps  were 
scaled  and  used  as  control.  A  total  of  31  horizontal  and  13  vertical  con¬ 
trol  points  were  used.  Standard  diviations  of  12  meters  in  X  and  Y  and  2 
meters  in  Z  (elevation)  were  assigned  to  the  control  for  weighting  purposes. 

The  74  points  in  the  model  were  photogrammetrically  aerotri- 
angulated  using  the  same  analytical  least  squares  solution  employed  for 
the  photography  flown  simultaneously  with  the  radar.  Data  reduction,  in¬ 
cluding  error  propagation,  was  performed  on  Autometric's  IBM  360/44  com¬ 
puter.  The  aerotriangulated  positions  of  the  points  in  the  test  area  are 
given  in  Table  II.  It  is  quite  apparent  from  fable  II  that  there  is  only 
modest  variation  in  elevation,  the  range  between  highest  and  lowest  points 
teing  90  meters. 

The  a  posteriori  variance  of  unit  weight,  computed  fron  the 
least  squares  solution,  indicated  a  standard  deviation  of  0.040  mm  in  the 
measurement  of  each  image  coordinate,  in  agreement  with  the  value  obtained 
from  the  aerotr iangulation  of  the  photography  flown  simultaneously  with 
the  radar.  The  error  propagation  yielded  standard  deviations  of  6.8  meters, 
6.7  meters  and  8.7  meters  in  the  triangulated  X,  Y  and  Z  coordinates, 
respectively.  The  fit  of  the  triangulated  model  to  the  control  scaled 
from  the  topographic  maps  resulted  in  RMS  (Root  Mean  Square)  discrepancies 
of  8.6  meters,  7.3  meters  and  7.4  meters  in  the  X,  Y  and  Z-coordinates , 
respectively. 
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TABLE  II.  AEROTRIANGULATED  POSITIONS  OF  POINTS  IN  TEST  AREA 


Point 

No. 

Type 

(See  Code 

In  Radar 
)  Model 

UTM  Coordinates 

Elevation  Above 

X  (meters) 

Y  (meters) 

■  Mean  Sea  Level 
Z  (meters) 

1 

A 

I  .  H 

72  7472.3 

37  30  )  1  3.7 

301.1 

? 

A 

I. II 

729561.3 

3731481.3 

307.4 

3 

A 

I  t  I  I 

729643.2 

3731 773.8 

2  74 .6 

4 

A 

I  f  I  I 

729426.8 

3731979.3 

290.4 

5 

A 

I  V  I  I 

727185.4 

37  32608.7 

2  9  2.8 

6 

A 

i  t  I  I 

728385. 1 

3733027.4 

300.5 

7 

A 

727570.7 

3732582.9 

?6o.6 

8 

A 

1 

726558.1 

I  3732783.9 

278.6 

1  l 

8 

I 

727280.7 

3733179.3 

254.2 

12 

* 

!  i 

727575.8 

373  331  7.8 

2  78.2 

1  3 

A 

1. 1 1 

7?  7920.8 

3733196.0 

2  71.7 

14 

A 

i 

728984.5 

3736560.7 

2  9  3  .7 

15 

A 

i 

728023.9 

37  36515.5 

251  .0 

16 

A 

i 

,  727096.9 

3736604.2 

245.3 

1  7 

A 

i 

727316.1 

3737168.7 

251.4 

18 

A 

i 

726915.9 

37  37  1  39. 8 

2  34.8 

1 9 

C 

i*  1 1 

729338.8 

3738427.5 

242 .4 

20 

A 

i 

72  m  7  06 . 0 

3738204.0 

242 .4 

21 

C 

i 

729205.0 

3738142.6 

24  1.5 

22 

C 

i 

727291.0 

37378)2.5 

2  2  3.4 

23 

i 

726573.4 

3737473,4 

21  7.2 

24 

A 

1. 1 1 

729166.5 

3733486.6 

298.7 

25 

A 

i .  1 1 

728565.2 

3733431.5 

2  8  o  .  8 

26 

t 

i .i  i 

728781.5 

37  34005.2 

254.6 

27 

A 

i 

7270  78.0 

3739014.1 

218.4 

28 

A 

i 

726309.6 

3738892.9 

2  36.3 

24 

A 

i 

726697.7 

37343)7.7 

2  2  3  . 1 

30 

A 

i 

727560.9 

3  7  39  72  7.1 

2  2  5.5 

n 

A 

i 

726934.2 

37  3  4  8  •>  3  .  4 

222.4 

32 

A 

i 

727244.7 

37401 39. 0 

2  3  7. 9 

i  3 

A 

i 

7263  3  1 .6 

3740824. 3 

302 . 7 

34 

i 

,  728034.5 

3740226. 7 

24  7.3 

36 

•- 

i,ii 

724126. 3 

3740745. 1 

2  3o.l 

36 

c 

i 

728468.9 

37  3  92  ,4.2 

21  ).6 

38 

G 

i .  1 1 

7  2  7451 .4 

3738559.2 

2  3  A  .  5 

39 

A 

i .  1 1 

72  774  7.4 

3740554. 3 

26  3  ,, 

4  3 

1- 

i 

729714.8 

3734158.3 

2  >7.0 

44 

A 

i .  1 1 

i  728485.5 

37  34  326.8 

2  5  8.9 

4  5 

A 

i .  1 1 

727803.9 

373,116.3 

25  )  .4 

46 

A 

i 

7?  7829.4 

3  7  354  70. 1 

?3o.3 

48 

[3 

1 1 

7  30662.0 

3734571.6 

2  4  .5 . 0 

49 

13 

i ,  n 

724  356. 8 

3736714.6 

239.6 

60 

B 

i .  1 1 

72  7594. 1 

37  36744.6 

2  34.2 

51 

A 

1 1 

72  7904. 5 

3736641,8 

25  7.0 

6  3 

A 

1 1 

7  301  88.2 

3737470.5 

254.4 

>7 

G 

i .  1 1 

728879. 3 

37402  70.6 

244.3 

58 

C 

•  i 

730103.1 

3  74060  7.5 

260.6 

59 

D 

1 1 

730682.9 

3740541.0 

248.1 

61 

A 

i .  1 1 

728650.  9 

3740803.5 

2  50.6 

62 

A 

i .  1 1 

728345.9 

37406  7  3.4 

265.6 

63 

(- 

i .  1 1 

728600.9 

3740199.3 

2  36.2 

64 

A 

i .  1 1 

7299  30.5 

3739970.5 

2  34.7 

90 

A 

1. 1 1 

729457.6 

3730973.6 

301.7 

9  6 

A 

i .  1 1 

72  7994.7 

- 

3734426.4 

2  7  7.0 

Code:  A  Road  Intersection  E  Point  of  Trees 

B  Point  of  Land  in  Lake  F  Wooded  Fence  Corner 

C  Building  Corner  G  Cupola  on  Building 

D  Intersection  of  Road  and  Drainage  Feature 
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Radar  Data  Processing  (Task  7) 

2*  ’ 1  Preliminary  Stereo  Radar  Test  Runs 

■  Hr  ^SlnC£  tbe  cJuality  of  the  test  material  was  less  than  orig- 
special^nnrir^f  ’  3nd  lack  °f  certain  data  required  adoption  of  a 
doubt\hatPr  ^  dat3  proressin8>  thare  was  considerable 

under  t£ese ^  P°Si^ionS  determined  by  stereo  radar  intersection 
associated6  *  WOUld  be  meaningful‘  With  this  uncertainty 

radar  te^r  W1  6  pr°ject’  it:  was  decided  to  make  preliminary  stereo 
reduction  n^Lf3  S°°a  f  P°SSible  in  order  to  verify  the  adopted  data 
expected  in  the  ^  ^  determine  the  order  of  accuracy  which  might  be 
resets  could  tT  h"  intfSected  test  poi"t  positions.  The  preliminary 

continue  wirh  the  project  "  3  b3SiS  ^  deCidlag  whether  °r  "0t  C° 

(oDDosite  lir  prel^min^y  runs  were  first  made  for  Stereo  Radar  Model  I 
terms  used  in^336?  ‘  ^  Varied  with  aspect  to  the  number  of 

to  the  nul  f6  5  ^  ^  polynomials.  and  with  respect 

tha^  rh  I  of  ground  points  used  for  control.  The  runs  revealed 
that  the  quadratic  term  could  be  dropped  from  the  slant  range  polynomial 

PO  nts  couLl^^d033;  1058  in  and  thG  number  of8contror 

points  could  be  reduced  considerably  with  very  little  loss  in  accuracy 

reduct ion'nroc 'd  °WeVer’  Pb*  encouraging  revelation  that  the  data"' 
reduction  procedure  worked  and  the  general  level  of  accuracy  of  the 

of  the^adar^ntf  s^t  d°int  P°Siti°nS  WaS  surPrisingly  good.  Comparison 
_  e  radar  intersected  test  point  positions  with  corresponding  aero- 

trisngulated  positions  resulted  in  RMS  (Root  Mean  Square)  discrepancies 

Inhr  lT;:dl,fT  7;6  t0.9-6  "eters  ln  X-  9'8  a>’l  -ter.  i"r 

and  11. U  to  19.3  meters  in  Z . 

A  sixth  preliminary  run  was  made  for  Stereo  Radar  Model  II 
(same-side  case),  with  resulting  RMS  of  discrepancies  for  the  test 
points  of  13.0  meters  in  X,  12.9  meters  in  Y,  and  20.1  meters  in  Z. 

on  January T30Sei97^ti3T\reSUlJS1WerK  r6vieWed  ^  Government  representatives 
,y  ’  1970 •  11  was  felt  that  the  results  were  sufficiently 
encouraging  to  warrant  continuation  of  the  project  to  completion. 

Additional  preliminary  runs  for  both  stereo  radar  models  were 
msde  to  check  the  accuracy  using  very  few  control  points,  and  to  edit 
out  any  bad  data.  The  accuracy  held  up  remarkably  well  when  very  few 
control  points  were  used,  and  only  four  points  exhibits  gross  dLcrep- 

M27  TdYnY  °f  th!lr  co°rdl"«es-  Three  of  the  bad  points 

Jell’inuLTl  n 5°,.  ^1>eared  ln  Model  U  two  of  the  points  (#58  and  1/90) 
lJ.  4°de1.11,  .  hese  points  were  rejected  from  further  consideration 
leaving.  4b  points  in  Model  I  and  28  poims  in  Model  II  for  control  and  ' 

runs.  ^  St3ge  ^  ^  S6t  f°r  or8anizin8  a  series  of  final  test 
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2.7.2 


Final  Stereo  Radar  Test  Runs 


The  final  stereo  radar  test  runs  were  organized  for  the 
purpose  of  revealing  the  effects,  if  any,  of  variation  in  the  reduction 
technique,  variation  in  control  point  configuration,  and  variation  in 
radar  image  coordinate  weighting. 

Variation  in  the  reduction  technique  was  made  with  respect 
to  definition  of  the  radar  flight  path,  evaluation  of  the  BgR-coefficient 
in  the  slant  range  calibration,  and  definition  of  the  zero-Doppler  line. 
The  radar  flight  path  for  each  strip  was  defined  either  by  a  set  of 
line  segments  connecting  successive  camera  stations,  as  given  in 
Section  2.3.1,  or  by  a  single  straight  line  between  the  first  and 
.Last  camera  stations.  The  slant  range  calibration  was  based  upon 
whether  the  BgR-coef ficient ,  appearing  in  Equation  (7)  (Section  2.3.3), 
was  determined  separately  for  each  record,  or  from  all  three  records 
combined.  (The  CgR-coef ficient  was  dropped  from  the  polynomial  in 
accordance  with  results  of  the  preliminary  test  runs  which  indicated 
that  the  quadratic  term  was  unnecessary) .  The  zero-Doppler  line  was 
defined  either  normal  to  the  flight  path,  as  it  should  be,  or  normal 
to  the  ground  track,  as  a  simplified  model  might  have  it. 

Four  distinct  control  point  configurations,  shown  in  Figure  12, 
were  employed  in  the  final  test  runs.  Basically,  14  points,  distributed 
fairly  evenly  over  the  test  area,  were  selected  as  the  control  points 
for  both  stereo  radar  models  (Model  I:  Opposite-side  case;  and  Model  II: 
Same-side  case).  In  Configuration  (A),  all  14  points  were  used  to 
calibrate  slant  range  and  correlate  "time".  in  Configuration  (B) ,  all 
14  points  were  again  used  for  the  slant  range  calibration,  but  only  2 
of  these  points,  one  near  each  end  of  the  test  area,  were  used  to  corre¬ 
late  "time".  In  Configuration  (C) ,  only  7  of  the  points  were  used  for 
the  slant  range  calibration  and  the  "time"  correlation.  Three  points 
were  selected  at  each  end  of  the  test  area,  and  the  7th  point  was  located 
near  the  middle  of  the  area.  In  Configuration  (D) ,  only  4  points,  two 
at  each  end  of  the  test  area,  were  used  to  calibrate  slant  range,  while 
2  of  these  points  (the  same  2  used  in  Configuration  (3))  were  used  to 
correlate  "time". 

For  the  most  part,  weighting  of  the  radar  image  coordinates, 

RI  and  Ti,  was  uniform;  that  is,  a  standard  deviation  of  0.030  mm  was 
assigned  to  each  Rj  and  to  each  Tj.  For  one  test  run,  however,  Rx 
was  heavily  weighted  (assigned  standard  deviation  =  0.0001  mm),  while 
the  weight  of  Tx  remained  unchanged  (0.050  mm);  for  another  test  run, 

Tx  was  heavily  weighted,  while  the  weight  of  Rx  remained  unchanged. 

The  series  of  final  stereo  radar  test  runs  is  comprised  of  11 
runs  for  each  of  Stereo  Radar  Models  I  and  II.  The  parameters  for  each 
of  the  11  runs  are  listed  in  Table  III. 

Since  14  points  were  selected  for  control,  32  points  remained 
in  Model  I  and  14  points  remained  in  Model  II  as  test  points. 

The  computer  programs  1  through  4,  described  in  Section  2.3.7, 
were  used  in  sequence  to  process  the  series  of  test  runs.  Then,  for 
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FIGURE  12  CONTROL  POINT  CONFIGURATIONS 


TABLE  III.  PARAMETERS  FOR  THE  FINAL  STEREO  RADAR  TEST  RUNS 


35 


r:;.h  tCtt  rmJ»  final  radar  intersected  positions  of  the  test  points 

ind*’rhe,nRM^-Wlth  the.corresP°nding  aerotriangulated  positions  (Table  II), 
.  discrepancies  were  computed  for  each  of  the  three  coordinates, 

‘  (  L r oss-track) ,  Y (along-track) ,  and  Z(elevation) .  The  RMS  values  are 
documented  later  in  Section  2.8. 

2,7,3  IiL  °f  Radar  Records  to  Ground  Points 


In  addition  to  determining  planimetric  position  and  elevation 
-rom  a  pair  of  radar  records  using  stereo  radar  techniques,  it  is  also 
possible  to  determine  planimetric  position,  but  not  elevation,  directly 
from  a  single  record.  This  approach  is,  of  course,  not  as  rigorous 
s  the  stereo  radar  approach,  since  it  does  not  account  for  the  effects 

f  dlsPdacement-  However,  in  view  of  the  modest  relief  actually 

exhibited  m  the  selected  test  area  (about  90  meters  maximum),  it  was 
felt  that  relief  displacement  would  not  be  very  detrimental. 


Essentially,  what  is  done  is  to  fit  the  radar  record  to  the 
aerotriangulated  positions  using  a  coordinate  transformation.  If 
the  radar  record  is  indeed  a  faithful  geometrical  representation,  in 
plan,  of  the  terrain  it  is  recording,  then  a  good  fit  should  result. 
Otherwise,  large  residuals  will  occur  at  the  fitted  points 


Since  the  radar  records  for  this  project  are  presented  in  slant 
range,  conversion  to  ground  range  presentation  is  necessary  before  a 
faithful  planimetric  representation  of  the  terrain  can  be  obtained. 

Ground  ranges  were  computed  using  the  following  equation: 


(16) 


In  Equation  (16) ,  GR  is  the  ground  range  presentation  at  image  scale; 

Kl  is  the  plate  coordinate,  as  computed  from  Equation  (5),  Section  2.3.3- 

u  t  fly1"8  ^18ht  above  terrain,  given  by  Lhe  aerotriangulation  of 
he  photography  flown  simultaneously  with  the  radar;  and  A™  and  Bcu 
are  parameters  of  the  slant  range  calibration  polynomial  (Equation  (7) 
m  Section  2.3.3). 


Although  the  slant  range  presentation  is  not  really  a  faithful 
planimetric  representation  of  the  terrain,  it  was  decided  to  include 
fitting  the  slant  range  presentation  to  the  ground  control  anyway  to 
see  how  much  it  affected  the  fit.  Hence,  both  ground  range  and  slant 
range  presentations  are  fitted. 


Two  different  coordinate  transformations  were  used  for  fitting 
the  data.  The  first  transformation  is  given  by  Equation  (17): 
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(17) 


lhis  is  a  4-parameter  transformation;  the  4  parameters  are  an  overall 
scale  factor,  X,  a  rotation,  (f>,  and  two  translation  constants,  Ki  and  K2 . 
Rj  and  are  the  plate  coordinates  of  the  radar  image  of  point  I,  as 
computed  from  Equations  (5)  and  (6) ,  respectively;  and  Xj  and  Yj  are 
the  across-track  and  along-track  ground  coordinates,  respectively. 


The  second  transformation  which  was  used  is  a  6-parameter 
transformation,  given  by  Equation  (18): 


cos  cj> 


LIJ 


-sin  (j> 


sin  cp 


cos  <p 


1  sin  n 
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(18)  . 


In  this  transformation,  two  distinct  scale  factors,  X  and  X  ,  are 
employed  for  scaling  Rj  and  Tj ,  respectively;  and  a  sixth  parameter,  r), 
accounts  for  non-orthogonality  of  the  Rj.  and  Tj  axes. 


All  ground  points  measured  in  each  radar  record  and  used  in 
the  stereo  radar  test  runs  were  used  for  fitting  the  record  to  the 
aerotriangulated  positions.  Specifically,  52  ground  points  were  used 
to  fit  Strip  3  to  the  aerotriangulated  positions,  46  ground  points  were 
used  to  fit  Strip  4,  and  28  ground  points  were  used  to  fit  Strip  5. 


Altogether,  four  runs  were  made  for  each  radar  strip.  Table  IV 
indicates  what  was  done  for  each  run. 


TABLE  IV.  FIT  OF  RADAR  RECORDS  TO  GROUND  POINTS 


Run 

No. 

Presentation  of 

Radar  Imagery 

Type  of 

Transformation 

1 

Uncorrected  Slant  Range 

4-Parameter 

2 

Uncorrected  Slant  Range 

6-Parameter 

3 

Corrected  to  Ground  Range 

4-Parameter 

4 

Corrected  to  Ground  Range 

6-Parameter 

In  each  case,  the  residuals  at  the  ground  points  were  deter¬ 
mined,  and  RMS  values  for  X  and  Y  were  computed  from  the  residuals.  The 
RMS  values  are  shown  in  Table  X,  Section  2.8.2. 

2.8  Analysis  of  Results  (Task  8) 

2.8.1  Topographic  Accuracy  of  Stereo  Radar 

The  results  of  the  final  stereo  radar  test  runs,  the  parameters 
of  which  are  indicated  in  Table  III,  are  given  in  Table  V.  These  results 
consist  of  RMS  values  in  X,  Y  and  Z  for  each  of  the  two  stereo  radar 
models,  computed  from  discrepancies  between  the  radar  determined  positions 
and  corresponding  aerotriangulated  positions.  For  Stereo  Radar  Model  1, 
the  RMS  values  are  based  upon  32  test  points;  for  Stereo  Radar  Model  II* 
the  RMS  values  are  based  upon  14  test  points. 
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TABLE  V.  RESULTS  OF  THE  FINAL  STEREO  RADAR  TEST  RUNS 


Note.  (1)  Average  RMS  values  are  computed  from  Runs  1-10. 

(2)  X  is  in  the  across-track  direction,  Y  is  in  the 
along-track  direction,  and  Z  is  elevation. 


Except  for  Run  11,  the  results  are  very  good  and  very  encour¬ 
aging.  particularly  when  a  considerable  portion  of  the  error  indicated 
v  these  RMS  values  might  be  attributed  to  the  standard  for  comparison 
(the  aero triangulated  positions  exhihit  RMS  errors  of  8.6  meters, 

.--  meters,  and  7.4  meters  in  X,  Y  and  Z,  respectively,  when  compared 
W^t  corresP°nc^n8  map-scaled  positions).  Average  values,  representative 
of  the  results  and  computed  from  all  test  runs  except  No.  11,  are  in¬ 
cluded  at  the  bottom  of  Table  V.  From  these  averages,  two  very  general 
remarks  can  be  made: 

(1)  Model  I  (Opposite-side  case)  provides  greater  accuracy 
than  Model  II  (Same-side  case),  particularly  in  Y, 

the  along-track  coordinate. 

(2)  The  X  (across-track)  coordinate  is  determined  with 
greater  accuracy  than  the  Y  and  Z  coordinates  for 
both  Models  I  and  II. 

The  discrepancies  in  planimetry  for  Run  4  are  plotted  in 
vector  form  in  Figures  13  (Model  I)  and  14  (Model  II).  These  figures 
show  how  the  discrepancies  are  distributed  over  the  test  area,  in 
planimetry  at  least.  It  is  apparent  from  these  plotted  discrepancies 
that  certain  systematic  effects  are  present.  An  attempt  was  made  to 
ascertain  the  cause  or  causes  of  the  systematic  effects,  but  it  was 
without  success. 

In  order  to  demonstrate  the  effects,  if  any,  of  processing 
t  e  data  under  different  conditions,  the  test  runs  are  re-grouped  into 
four  specific  sequences.  The  first  sequence  (Runs  1,  2,  3  and  4)  at¬ 
tempts  to  show  the  effects  of  varying  the  reduction  technique  in  two 
ways : 

(1)  The  flight  path  is  defined  by  a  set  of  line  segments 
connecting  successive  camera  air  stations,  as 
described  in  Section  2.3.1;  or  it  is  defined  as  one 
straight,  line  between  the  first  and  last  stations. 

(2)  The  BSR-coefficient  of  the  slant  range  calibration 
polynomial  (Equation  (7))  is  evaluated  separately 
for  each  radar  pass  (3,  4  or  5);  or  it  is  evaluated 
only  once  using  data  from  all  3  passes  combined. 

For  convenience  in  making  comparisons,  the  results  of  this  sequence 
of  runs  are  repeated  in  Table  VI.  All  parameters  other  than  those  under 
investigation  are  held  fixed  (See  "Note"  under  Table  VI). 

The  second  sequence  (Runs  4,  5,  8  and  9)  is  designed  to  show 
the  effect  of  varying  the  reduction  fc  'nique  with  respect  to  how  the 
zero-Doppler  line  is  defined;  that  if-  normal  to  the  flight  path,  or 
normal  to  the  ground  track.  The  results  of  this  sequence  are  repeated 
in  Table  VII.  Tests  were  made  using  two  different  ground  control 
configurations,  A  and  D,  Figure  12. 
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FIGURE  14  DISCREPANCY  VECTORS  IN  PLANIMETRY 
FOR  STEREO  RADAR  MODELS  (RUN  4). 


TABLE  VI.  TOPOGRAPHIC  ACCURACY  OF  STEREO  RADAR:  EFFECT  OF 

FLIGHT  PATH  DEFINITION  AND  B  -COEFFICIENT  EVALUATION 

oK 


RMS  VALUES 

Run 

Flight 

Evaluation 

Model  I 

Opposite-side 

Model 

II:  Same-side 

No. 

Path 

...  of  bsr 

X(m) 

Y(m) 

Z(m) 

X(m) 

Y(m) 

Z(m) 

1 

Segmented 

Separate 

11.9 

10.0 

14.7 

10.6 

14.7 

14.1 

2 

Segmented 

Combined 

9.5 

9.9 

11.0 

11.3 

14.7 

11.8 

3 

Straight 

Line 

Separate 

7.4 

11.8 

13.9 

9.0 

21.4 

16.2 

4 

Straight 

Line 

Combined 

7.3 

11.8 

12.3 

8.9 

21.5 

15.4 

Note : 


(1)  All  runs  in  this  table  use  ground  control 
configuration  A  (14  points  for  slant  range 
calibration  and  for  time  correlation), 

(2)  The  zero-D oppler  line  is  defined  normal  to 
the  flight  path. 

(3)  Rj  and  T^  are  equally  weighted. 
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TABLE  VII.  TOPOGRAPHIC  ACCURACY  OF  STEREO  RADAR: 

EFFECT  OF  ZERO-DOPPLER  LINE  DEFINITION 


Run 

No. 

Zero-Doppler 
Line  Normal 

To 

Ground 

Cont  rol 
Configuration 

RMS  Values 

Model 

I : Opposite-side 

Model 

II : Same-side 

X(m) 

Y  (m) 

Z(m) 

X(m) 

Y(m) 

Z(m) 

r  4 

Flight  Path 

A 

7.3 

11.8 

12.3 

8.9 

21.5 

15 . 4 

5 

Ground  Track 

A 

7.3 

11.8 

12.3 

8.7 

21.5 

13.5 

8 

Flight  Path 

D 

6 . 5 

1 3  3 

14.8 

10.6 

21,1 

24.8 

9 

Ground  Track 

D 

6.5 

13.3 

14.9 

9.9 

21.1 

22.5 

Note.  (1)  Ground  control  configuration  A  uses  14  points 

for  slant  range  calibration  and  time  correlation; 
configuration  D  uses  4  points  for  slant  range 
calibration  and  2  for  time  correlation. 

(2)  For  all  runs  in  this  table,  the  flight  path 
is  defined  as  a  straight  line,  and  BSR  is 
evaluated  from  all  j  records  combined. 

(3)  Rj  and  Tj.  are  equally  weighted. 
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The  third  sequence  (Runs  4,  6,  7  and  8),  repeated  in  Table  VIII, 
s  lows  the  effect  of  using  the  different  ground  control  configurations 

A,  B,  C  and  D(See  Figure  12),  while  holding  all  other  parameters  con¬ 
stant  . 


The  fourth,  and  final,  sequence  (Runs  4,  10  and  11),  repeated 
in  Table  XX,  shows  the  effect  of  image  coordinate  weighting.  Again, 
all  parameters  other  than  those  under  investigation  are  held  fixed. 

It  is  quite  evident  from  the  RMS  values  listed  in  Tables  VI, 
VII  and  VIII  that  no  drastic  changes  in  accuracy  take  place  with 
respect  to  the  variations  made  in  the  reduction  technique  anu  in  the 
ground  control  configuration.  This,  perhaps,  should  not  be  too  sur¬ 
prising  so  long  as  the  aircraft  flies  a  straight  and  level  course,  the 
radar  equipment  remains  stable  in  its  image  presentation,  and  the  test 
area  is  relatively  small  conditions  which  are  either  true  (the  test 
area  is  indeed  small)  or  very  probable. 


Nevertheless,  there  are  indications  in  the  tables  that  the 
accuracy  responds,  in  part,  to  some  of  the  changes.  For  example, 

Table  VI  shows  that  segmentation  of  the  flight  path  tends  to  decrease 
the  accuracy  of  the  X-coordinate  while  improving  the  accuracy  of  the 
Y-coordinate ,  particularly  in  Model  II;  and  that  the  accuracy  of  the 
Z-coordinate  tends  to  improve  when  all  3  radar  records  are  combined 
to  evaluate  the  Bg^-coef f icient .  Also,  Table  VIII  shows  that  the 
accuracy  of  the  Z-coordinate  in  Model  II  decreases  as  fewer  ground 
control  points  are  used  for  slant  range  calibration. 

Notwithstanding  the  relative  smallness  of  the  test  area  used 
in  this  sfudy ,  it  is  still  heartening  to  observe  in  Table  VIII  that  the 
few  ground  points  used  in  Configuration  D  do  about  as  good  a  job  for 
Radar  Model  I.  as  all  14  points  used  in  Configuration  A,  and  it  is  only 
in  the  Z-coordinate  of  Model  II  that  any  deterioration  in  accuracy  takes 
place.  Indeed,  it  seems  remarkable  that  such  accuracy  as  indicated  for 
Run  8  has  been  achieved  by  using  only  4  control  points,  particularly 
when  much  of  the  error  could  be  attributed  to  the  aerotriangu] ated 
positions  used  as  the  standard  for  comparison. 

ut  is  quite  obvious  from  Table  IX  that.  ueavily  weighting  the 
Rj  image  coordinate  has  very  little  erfect,  if  any,  on  accuracy,  perhaps 
a  slight  improvement  in  the  Z-coordinate  of  Model  II.  On  the  other 
hand,  heavily  weighting  the  Tj  imap  ^ordinate  is  disastrous.  This 
severe  loss  of  accuracy  results  rr.  -nforcing  the  intersection  of  the 
two  zero-Doppler  lines  (which  lie  in  uear-parallel  planes)  at  the 
expense  of  the  relatively  well-determined  slant  ranges.  Actually, 
there  is  no  good  reason  in  evidence  at  this  time  to  weight  the  R,  and 
TI  image  coordinates  unequally.  Hence,  the  Run  11  results,  although 
interesting,  are  of  little  consequence. 

2.8.2  Geomettic  Analysis 

The  results  of  fitting  the  separate  radar  records  to  the  aero- 
triangulated  positions,  as  descrioed  in  Section  2,7.3,  are  given  in 
Table  X.  The  first  portion  of  this  table  repeats  the  information  already 
given  in  Table  IV. 
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TABLE  VIII.  TOPOGRAPHIC  ACCURACY  OF  STEREO  RADAR:  EFFECT  OF 
VARIOUS  GROUND  CONTROL  CONFIGURATIONS 
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TABLE  IX.  TOPOGRAPHIC  ACCURACY  OF  STEREO  RADAR: 

EFFECT  OF  IMAGE  COORDINATE  WEIGHTING 


Weights  of  Image 


RMS  Values 


^ooramat  es 

Rj  and 

Model 

Opposite-side 

Model 

II:  Same-side 

X(m) 

Y(m) 

Z(m) 

X(m) 

Y  (m) 

Z(m) 

Equal 

7.3 

12.3 

8.9 

21.5 

15.4 

Rj  heavily  weighted 

7.3 

12.3 

8.  7 

21.4 

13.4 

Tj  heavily  we  ft  d 

663. 6 

2684.7 

1302.1 

50.2 

3449.9 

For  all  runs  in  this  table,  the  flight  path 
is  defined  as  a  straight  line,  the  B  -coefficient 
is  evaluated  from  all  3  records  combined,  the  zero 
Doppler  line  is  defined  normal  to  the  flight  path, 
and  ground  control  configuration  A  is  used. 
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TABT.E  X.  RESULTS  OF  FITTING  RADAR  RECORDS  TO  GROUND  POINTS 


Type 

of 

Transformation 

RMS  Values 

Run 

No. 

Presentation  of 
Radar  Imagery 

Strip  3 

Strip 

4 

Strip  5 

X(m) 

Y(m) 

X(m) 

Y(m) 

X(m) 

Y(m) 

1 

Uncorrected  Slant 
Range 

4-parameter 

198.2 

59.7 

118.6 

47.7 

27.8 

34.3 

2 

Uncorrected  Slant 
Range 

6-parameter 

24.9 

12.4 

15.4 

20.0 

10.0 

29.4 

3 

Corrected  to 

Ground  Range 

4-parameter 

17.7 

14.3 

14.8 

20.3 

11.3 

34.6 

4 

Corrected  to 

Ground  Range 

6-parameter 

16.1 

12.4 

13.5 

20.0 

10.1 

29.4 
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As  might  be  expected,  Table  X  indicates  a  very  poor  fit  of  the 
radar  imagery  to  the  aerotriangulated  positions  when  uncorrected  slant 
ranges  and  a  4-parameter  transformation  are  used,  particularly  in  the 
\  (across-track,  or  slant  range)  direction  for  Strips  3  and  4.  For  Strip 
5,  the  fit  in  X  is  much  better,  reflecting  the  closer  agreement  between 
slant  range  and  ground  range  at  the  greater  offset  distance  (18.4  NM, 
as  compared  to  6.6  NM  and  9.2  NM  for  Strips  3  and  4,  respectively). 

Table  X  also  shows  that  the  fit  is  tremendously  improved  when 
a  6-parameter  transformation  is  employed,  even  when  slant  ranges  remain 
uncorrected.  This  improvement  takes  place  because  the  6-parameter  trans¬ 
formation  can  accommodate,  in  an  average  sense,  the  difference  in  scale 
in  X  between  slant  range  and  ground  range,  while  at  the  same  time  inde¬ 
pendently  fitting  to  the  "time"  scale  in  Y.  The  ability  of  the  6-parameter 
transformation  to  account  for  non-orthogonality  of  coordinate  axes  may 
also  contribute  to  the  improvement. 

For  Strip  5,  use  of  a  6-parameter  transformation  with  uncorrected 
slant  ranges  results  in  a  fit  every  bit  as  good  as  the  one  using  corrected 
ground  ranges  ,  reflecting  the  very  small  change  in  slant  range  scale  with 
respect  to  ground  range  at  the  large  (18.4  NM)  offset  distance. 

Quite  naturally,  Table  X  indicates  considerable  improvement  in 
the  fit  after  slant  range  measurements  are  corrected  to  ground  ranges. 

This  is  as  it  should  be,  because  the  fit  is  made  to  ground  control 
points  which  have  essentially  a  "ground  range  presentation".  As  might 
be  expected,  correction  to  ground  range  combined  with  a  6-parameter  trans¬ 
formation  gives  the  best  fit  of  all. 

Understandably,  the  RMS  values  in  X  obtained  from  geometric 
fits  of  the  separate  radar  records  to  ground  control  are  not  as  small 
as  their  stereo  radar  model  counterparts,  although,  quite  honestly, 
they  are  not  that  much  worse,  particularly  after  correction  is  first 
made  to  ground  range.  As  mentioned  before,  the  larger  RMS  values  result 
from  the  effects  of  relief  displacement.  Since  only  modest  relief  is 
present  in  the  test  area,  the  effects  are  very  small. 

It  is  interesting  to  observe  in  Runs  2,  3  and  4  of  Table  X 
that  although  the  RMS  in  Y  increases  with  increasing  offset  distance 
(offset  distances,  as  already  indicated  in  Table  I,  are  6.6  NM,  9.2  NM, 
and  18.4  NM  for  Strips  3,  4  and  5,  respectively)  the  RMS  in  X  actually 
decreases.  This  seems  to  reflect  the  decreasing  influence  of  relief 
displacement  on  X  as  the  offset  distance  increases,  while  at  the  same 
time  suggesting  that  the  acauracy  of  slant  range  measurement  is  in¬ 
dependent  of  offset  distance. 

The  discrepancies  at  the  individual  points  are  plotted  in 
vector  form  in  Figures  15,  16  and  17  for  Strips  3,  4  and  5,  respectively. 

It  is  clear  from  these  figures  that  certain  systematic  effects  are 
exerting  influence.  However,  as  with  the  discrepancies  plotted  in 
Figures  13  and  14,  attempts  to  determine  the  source  or  sources  of  the 
systematic  effects  were  unsuccessful. 
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Northing  in  Meters 


Northing  in  Meters 


FIGURE  17.  DISCREPANCY  VECTORS  IN  PLANIMETRY 
IMAGERY  FIT  FOR  RECORD  5. 


2 .  S .  3 


Slant  Range  Analysis 


The  primary  function  of  side— looking  radar  is  to  measure  slant 
range.  It  is  of  interest,  therefore,  to  check  how  well  this  function 
is  being  performed. 

Slant  range  calibration  curves  for  a  representative  test  run, 
Run  4,  are  plotted  in  Figure  18.  For  this  run,  14  control  points  were 
used  for  the  slant  range  calibration.  Similar,  almost  identical, 
graphs  can  also  be  plotted  for  the  other  test  runs.  Figure  18  clearly 
shows  the  linear  relationship  between  the  image  coordinate,  F.j ,  and  the 
slant  range,  indicating  a  geometrically  undistorted  presentation  of  slant 
range  across  the  radar  record,  and  justifying  omission  of  the  quadratic 
term  in  the  calibration  polynomial.  Figure  18  also  shows  the  constancy 
of  slope  from  one  record  to  the  next,  indicating  stability  of  the  radar 
in  recording  slant  range,  and  justifying  computation  of  one  slope  term 
(the  BgR-coef f icient)  for  all  three  records. 

The  small,  almost  imperceptible,  departures  of  the  plotted 
points  from  the  straight  lines  fitted  through  them  demonstrate  the  high 
degree  of  accuracy  which  can  be  achieved  in  the  calibration.  RMS  values 

computed  from  the  departures,  or  residuals,  obtained  for  Test  Run  4 
are : 


Strip  3:  6.4  meters 

Strip  4:  13.0  meters 

Strip  5:  5.9  meters. 

These  values  are  based  upon  all  14  control  points.  The  other  test  runs 
using  the  same  14  control  points  for  slant  range  calibration  yield 
essentially  the  same  results.  Test  runs  using  fewer  control  points  ex¬ 
hibit  much  smaller  RMS  values.  However,  these  smaller  values  probably 
resulted  from  the  fewer  degrees  of  freedom  in  fitting  the  curves,  and 
so  are  not  given  as  much  credence  as  the  values  listed  above. 


The  RMS  values  for  Strips  3  and  5  seem  to  agree  with  each  other, 
and  are  quite  reasonable  in  light  of  the  accuracies  obtained  for  the 
radar  determined  positions,  Table  V.  These  two  nearly  equal  values 
also  suggest  that  slant  range  accuracy  is  independent  of  distance,  since 
the  offset  distances  for  Strips  3  and  5  are  6.6  NM  and  18.4  NM,  respec¬ 
tively.  On  the  other  hand,  tne  RMS  value  for  Strip  4  seems  excessive, 
and  without  any  apparent  reason.  Certainly,  radar  image  quality  cannot 
be  blamed,  for,  in  this  respect,  Strip  4  is  best. 

2.8.4  Error  Propagation  Results 

The  error  propagation  procedure,  described  Vriefly  in  Section 
2.3.6,  was  used  to  compute  standard  deviations  in  X,  and  Z  for  each 
test  point  position  determined  by  stereo  radar  intersection.  The  error 
propagation  was  based  upon  a  standard  deviation  of  0.050  mm  assigned 
to  each  radar  image  coordinate.  The  standard  deviations  were  averaged 
so  as  to  provide  measures  of  precision  which  would  be  typical  of  the 
entire  test  area.  These  averaged  values  are  given  in  Table  XI. 
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TABLE  XI.  STEREO  RADAR  TEST  POINT  PRECISION 


Stereo  Radar  Model 

Averaged  Standard  Deviations 

ax(m) 

aY(m) 

'z(m) 

Model  I:  Opposite-side  case 

6.3 

5.3 

9.2 

Model  II:  Same-side  case 

13.6 

5.3 

26.5 

* Applicability  of  Stereo  Radar  to  Topographic  Mapping 

Of  prime  concern  in  this  project  is  the  applicability  of  the 
AN/APQ-102 (XA-2)  side-looking  radar  to  1:50,000  and  1:250,000  scale  topo¬ 
graphic  mapping,  particularly  from  the  standpoint  of  accuracy.  The 
standards  of  accuracy  for  military  maps  are  given  in  AMS  Technical 
Manual  S-l,  "Specifications  for  Military  Maps".  The  standards  covering 
1:50,000  and  1:250,000  scale  topographic  maps  are  repeated  herein  for 
convenience . 

(1)  Class  A  Maps 

(a)  Horizontal  accuracy:  90%  of  all  planimetric  features,  except 
those  unavoidably  displaced  by  exaggerated  size  of  symbols, 
are  located  within  0.02  inch  of  their  geographical  position 
as  referred  to  the  map  projection. 

(b)  Vertical  accuracy:  90%  of  all  contours  and  of  all  elevations 
of  points  interpolated  from  contours  are  accurate  to  within 
half  the  basic  contour  interval. 

(2)  Class  B  Maps 

Maps  that  fail  to  meet  the  standards  of  accuracy  for  Class  A,  but 

meet  the  following  standards: 

(a)  Horizontal  accuracy:  90%  of  all  planimetric  features,  except 
those  unavoidably  displaced  by  exaggerated  size  of  symbols,  are 
located  within  0.04  inch  of  their  geographic  position  as  referred 
to  the  map  projection. 

(b)  Vertical  accuracy:  90%  of  all  contours  and  of  all  elevations 
of  points  interpolated  from  contours  are  accurate  to  within  one 
contour  interval. 

(3)  Class  C-l  Maps 

Maps  that  fail  to  meet  the  standards  of  accuracy  for  Class  B,  but 

meet  the  following  standards: 

(a)  Horizontal  accuracy:  90%  of  all  planimetric  features  shown, 
except  those  unavoidably  displaced  by  exaggerated  size  of  symbols, 
are  located  within  0.08  inch  of  their  geographical  position  as 
referred  to  the  map  projection. 

(b)  Vertical  Accuracy:  90%  of  all  contours  and  of  all  elevations 
of  points  interpolated  from  contours  are  accurate  to  within  two 
contour  intervals. 

(4)  Class  C-2  Maps 

Maps  that  fail  to  meet  the  standards  of  accuracy  for  Class  C-l,  but 
must  be  retained  to  fill  possible  emergency  needs,  pending  replace¬ 
ment. 
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Applying  these  standards  specifically  to  1:50,000  and  1:250,000 
scale  topographic  maps  yields  the  accuracy  requirements,  in  meters,  listed 
in  Tables  XII  and  XIII. 


M  Three  stereo  radar  test  runs,  yielding  overall  "best",  "average 
and  worst"  accuracies,  respectively,  were  selected  from  Table  V  for  com¬ 
parison  with  the  standards  of  accuracy.  Test  Run  2,  the  "best",  uses  a 
segmented  flight  path,  14  control  points  for  slant  range  calibration  and 
time  correlation,  and  a  zero-Doppler  line  defined  normal  to  the  flight 
path.  Test  Run  4,  the  "average",  uses  a  straight  line  flight  path; 
otherwise,  it  is  the  same  as  Test  Run  2.  Test  Run  9,  the  "worst"  (Test 
Run  11  was  excluded  from  consideration),  uses  a  straight  line  flight 
path,  only  4  control  points  for  slant  range  calibration  and  2  for  time 
correlation,  and  a  zero-Doppler  line  normal  to  the  ground  track. 


The  planimetric  errors  of  the  test  points  were  computed  from 
their  X  and  Y  components.  The  percentages  of  points  lying  within  the 
standards  of  horizontal  accuracy  for  Class  A  and  Class  B  1:50,000  scale 
maps.  Table  XII,  were  computed  for  all  3  test  runs,  and  for  both  stereo 
radar  models.  These  percentages  are  listed  in  Table  XIV. 

The  percentages  of  points  lying  within  the  standards  of  vertical 
accuracy  for  Class  A,  B  and  C-l  1:50,000  and  .:250,000  scale  maps, 

Tables  XII  and  XIII,  were  also  computed.  These  percentages  are  also 
listed  in  Table  XIV. 

In  each  instance,  the  first  percentage  equal  to  or  exceeding 
90%  is  underlined  in  Table  XIV.  From  these  data,  the  following  state¬ 
ments  are  made : 

(1)  Stereo  Radar  Model  I  yields  sufficient  horizontal  accuracy  for 
Class  A  mapping  at  1:50,000  scale. 

(2)  Stereo  Radar  Model  II  yields  sufficient  horizontal  accuracy  for 
Class  B  mapping  at  1:50,000  scale,  or  for  Class  A  mapping  at 
1:250,000  scale. 

(3)  Stereo  Radar  Model  I  yields  sufficient  vertical  accuracy  for 
Class  B  mapping  at  1:50,000  scale  with  20-meter  contours,  or 
for  Class  C-l  mapping  at  1:50,000  scale  with  10-meter  contours, 
provided  ample  control  is  available. 

(4)  Stereo  Radar  Model  I  yields  sufficient  vertical  accuracy  for 
Class  A  mapping  at  1:250,000  scale  with  50-meter  contours,  or 
for  Class  B  mapping  at  1:250,000  scale  with  25-meter  contours. 

(5)  Stereo  Radar  Model  II  yields  sufficient  vertical  accuracy  only 
for  Class  C-l  mapping  at  1:50,000  scale  with  20-meter  contours. 

(6)  Stereo  Radar  Model  II  yields  sufficient  vertical  accuracy  for 
Class  A  mapping  at  1:250,000  scale  with  50-meter  contours,  or 
for  Class  B  mapping  at  1:250,000  scale  with  25-meter  contours, 

only  under  the  "best"  conditions  of  Test  Run  2.  Otherwise,  Model  II 
can  be  used  only  for  Class  A  mapping  with  100-meter  contours,  or 
for  Class  B  mapping  with  50-meter  contours,  or  for  Class  C-l  mapping 
with  25-meter  contours. 
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TABLE  XII.  SPECIFIC  STANDARDS  OF  ACCURACY  FOR  1:50  000 
SCALE  TOPOGRAPHIC  MAPS 


Horizontal  Accuracy 
(90%  Standard) 

Vertical  Accuracy  (90%  Standard) 
Contour  Interval 

5  meters 

10  meters 

20  meters 

Class  A  Maps 

25.4  m . 

2.5  m . 

5.0  m. 

10.0  m. 

Class  B  Maps 

50.8  m. 

5.0  m . 

10.0  m. 

20.0  m. 

Class  C-l  Maps 

101.6  m. 

10 .0  m. 

20.0  m. 

40.0  m. 

TABLE  XIII.  SPECIFIC  STANDARDS  OF  ACCURACY  FOR  1:250  000 
SCALE  TOPOGRAPHIC  MAPS 


Class  A  Maps 
Class  B  Maps 
Class  C-l  Maps 


Horizontal  Accuracy 
(90%  Standard) 

127  m. 

254  m. 

508  m. 


~Vert ical  Accuracy- (90%  Standard)' 
Contour  Interval 

25  meters  50  meters  xoo  meters 
12.5  m.  25  m.  50  m 


12.5  m . 
25  m . 

=;n  m 


50  m. 


50  m 
100  m. 
2"~ 


TABLE  XIV.  PERCENTAGE  OF  POINTS  WITHIN  SPECIFIED  LIMITS 

FOR  THREE  SELECTED  TEST  RUNS 


5 


There  is  little  doubt  that  the  stereo  radar,  particularly  Model  I, 
provides  sufficient  horizontal  accuracy  for  1:50,000  and  1:250,000  scale 
mapping.  It  is  the  vertical  accuracy  which  is  apparently  insufficient  in 
certain  cases,  and  is  therefore  the  controlling  factor.  In  order  to  por¬ 
tray  more  graphically  what  has  been  said  in  statements  (3)  through  (6) 
with  respect  to  vertical  accuracy,  Tables  XV  and  XVI  have  been  prepared. 

Si_at"ements  (3)  through  (6),  and  the  accompanying  information 
in  Tables  XV  and  XVI,  should  be  accepted  with  caution,  bearing  in  mind 
the  following  qualifications: 

(1)  Th ■=*  vertical  accuracies  obtained  in  the  test  runs 
are  with  respect  to  well— identified  points,  measured 
individually;  that  is,  the  points  are  not  interpola¬ 
ted  from  map  contours,  as  the  mapping  specifications 
prescribe. 

(2)  Comparison  was  made  with  aerotriangulated  positions 
which  are  not  entirely  error-free  themselves. 

Part  of  the  error  assigned  to  the  radar  determined 
positions  can  thus  be  attributed  to  error  in  the 
accuracy  standard  against  which  comparison  was  made. 

(3)  The  test  area  has  only  modest  relief;  that  is,  the 
testing  was  done  only  with  moderately  hilly  terrain, 
certainly  not  with  mountainous  terrain. 

The  preceding  qualifications  notwithstanding,  the  results  which 
have  been  obtained  in  this  study  are  very  heartening,  for  they  clearly 
demonstrate  a  degree  of  geometrical  fidelity  in  the  AN/APQ-102 (XA-2) 
radar  that  offers  much  promise  for  application  to  1:50,000  and  1:250  000 
scale  topographic  mapping. 
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TABLE  XV.  SUITABILITY  OF  STEREO  RADAR  FOR  1:50,000 
SCALE  TOPOGRAPHIC  MAPPING 


Contour  Interval 

5  meters 

10  meters 

20  meters 

Class  A  Maps 

- 

- 

- 

Class  B  Maps 

- 

- 

I* 

Class  C-l  Maps 

- 

I* 

I,  II 

I*:  Stereo  Radar  Model  I  is  suitable,  only  when  ample 

control  is  available. 

I, II:  Stereo  Radar  Models  I  and  II  are  suitable. 


TABLE  XVI.  SUITABILITY  OF  STEREO  RADAR  FOR  1:250,000 
SCALE  TOPOGRAPHIC  MAPPING 


Contour  Interval  ! 

25  meters 

50  meters 

100  meters 

Class  A  Maps 

- 

I, II* 

I, II 

Class  B  Maps 

I, II* 

I, II 

I, II 

Class  C-l  Maps 

I, II 

I, II 

I, II 

I, II*:  Stereo  Radar  Model  I  is  suitable;  Model  II  is 
suitable  only  under  "best"  conditions. 

I, II:  Stereo  Radar  Models  I  and  II  are  suitable. 
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DISCUSSION 


3 • 1  Test  Material 

No  particular  problems  were  encountered  in  sorting  out  the  test 
material.  The  lack  of  any  kind  of  indexing  for  the  vertical  photography 
did,  however,  make  the  job  of  identifying  the  imagery  by  strip  a  tedious 
one.  During  the  inspection,  some,  of  the  data  were  disqualified  from 
consideration  as  test  material: 

a)  Vertical  photography  of  Strips  1  and  2  for  both 
Sortie  08—61  and  Sortie  08-75  was  exposed  at 
60^  end lap .  Strips  having  more  closely  spaced 
camera  stations  (more  exposures)  can  better  de¬ 
fine  the  flight  path  for  the  purpose  of  the  in¬ 
vestigation. 

b)  Radar  Pass  6  of  Sortie  08—61  had  no  simultaneous 
vertical  photography. 

c)  Pass  4  of  Sortie  08-75  had  no  radar  imagery. 

d)  Image  quality  on  four  radar  strips  (08-61,  Pass  6; 

08-75,  Passes  3  and  5;  and  08-212,  Pass  3)  was 
very  poor. 

In  view  of  the  poor  quality  of  the  photography  and  radar  imagery 
obtained  for  this  investigation,  it  was  decided  by  Autometric,  with  the 
concurrence  of  ETL,  to  confine  the  test  to  an  area  smaller  than  that 
originally  planned.  Radar  image  quality  had  the  largest  influence  on  the 
selection  of  test  material.  Unfortunately,  the  vertical  photography 
corresponding  to  the  best  radar  coverage  was  of  poor  quality.  Radar 
image  quality  also  influenced  the  location  of  the  test  area,  since  the 
size  of  the  test  area  was  limited  and  image  quality  varied  within  the 
strip. 


3 . 2  Stereoscopic  Viewing 

During  the  course  of  the  data  inspection,  it  became  doubtful  to 
the  investigators  that  the  conjugate  radar  imagery  provided  could  be 
visually  correlated  to  the  point  where  stereoscopic  parallaxes  could  be 
measured  directly.  Only  slight  relief  was  apparent  in  the  same-side 
radar  model,  and  very  little  relief,  real  or  imagined,  was  seen  in  the 
opposite-side  stereomodel.  In  fact,  for  the  opposite-side  model,  a  good 
deal  of  visual  effort  was  required  just  to  correlate  images.  When  a 
floating  mark  was  introduced  into  the  model,  it  was  barely  possible,  if 
at  all,  to  maintain  correlation  of  imagery  and  perceive  the  floating  mark 
simultaneously.  This  is  understandable  since  the  principle  of  projection 
in  the  human  eye  is  different  from  that  of  the  radar  projection.  The 
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absence  of  noticeable  relief,  even  in  the  same-side  model  where  images 
were  easily  correlated,  can  then  be  explained  by  the  difference  between 
a  and  radar  parallax.  The  same-side  model  can  be  visually  correlated 
Decause  the  radar  parallaxes  are  much  smaller  than  those  of  the  opposite- 

However^hpra  the  Gyes '  ^ability  to  fuse  conjugate  images. 

*  becfUS?  of  the  smaHer  parallaxes,  little  relief  is  seen.  In 

ne^hair%a?^al  ground  reHef  in  the  area  is  only  200  to  300  feet, 
perhaps  yielding  radar  parallaxes  too  small  to  have  visual  significance. 


3.3 


Data  Reduction  Approach 


v-i.nal  -  nce  5  Was  0bvi0us  that  little  could  be  accomplished  using 
IntirPlv  re°SC?Pi  •  mfh0ds’  the  aPProach  in  the  investigation  relied 
" ly  °n  analytical  methods,  along  the  lines  of  analytical  photogram- 

The  mlthemathr°?8h  ^  ?e°mebry  and  conditions  imposed  by  the  radar  system. 

incor™^  ?  ySuS  Programming  for  the  investigation  had  to 

lb?  techniques  for  calibrating  the  radar  range  measure- 

a  ol  1  1118  311  a[tificial  time  scale ,  and  correlating  time,  since 

SL  th  I'1  data’  such  as  range  marks  and  time  marks,  were  missing 

dlta  in  s.ch  reCOrda‘  /S  “  tUrned  °Ut*  the  establishment  of  missing 
data  in  such  an  empirical  manner  was  quite  successful. 

±. j  .  In  mathematical  analysis,  the  original  thought  was  to  fit 

iiterno^pmhpr  ^  ^  fli8ht  P*th  (^^ngulated  camera  stations)  and 
interpolate  between  camera  stations  for  radar  positions  using  second  or 

higher  order  polynomials.  However,  with  the  camera  station  spacing  so 

Zfl  T  «-  larfiest  spacing  represlnts 

about  1/2  minute  of  flying  time,  it  seemed  impractical  to  fit  hieh-order 

polynomials  to  the  flight  path.  As  a  result,  linear  functions  were  used 

•  the  £llSht  paths,  with  2  successive  camera  stations  defin- 

mg  each  increment  of  the  flight. 


3.4 


Radar  Mensuration 


It;  waa  originally  intended  to  work  with  enlarged  scale  radar 
magery  (approximately,  four  times  negative  scale).  However  after  a 
visit  to  the  Engineering  Topographic  Laboratories  for  the  purpose  of 
viewing  the  f , larged  imagery  on  the  Wild  STK-1  Stereocomparator  it  be¬ 
came  evident  that  the  additional  enlargement  from  optica?  magnification 
in  the  instrument  degraded  the  model.  Since  the  Mann  Comparator  has 
simlar  optical  magnification,  it  was  decided  to  use  imagery  at  negative 


Mf.  ®  ma?or  dlfflculty  in  the  radar  mensuration  task  was  the  id 

tificrition.of  points  on  the  radar  imagery.  Image  quality  was  lowest  in"' 
the  same-side  model  where  only  30  points  were  identified  for  measurement 
The  opposite-side  model  provided  an  adequate  total  of  49  points 


3.5 


Aerotriangulation 


Save^al  Problems  were  encountered  in  the  aerotriangulation  of 
h  vertical  photography  taken  simultaneously  with  the  radar  The  1  i 
of  camera  calibration  parameter  made  it  impossible tl  compensate  fo^ 
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lilra  distortion,  and  the  minimal  fiducial  imagery  hampered  the  transfor¬ 
mation  of  comparator  measurements  to  plate  coordinates.  It  was  necessary 
to  define  the  plate  coordinate  system  and  assume  the  location  of  the 
principal  point  from  comparator  measurements  of  the  three  available 
fiducial  marks.  The  principal  point  was  determined  by  computing  the  in¬ 
tersection  of  the  line  formed  by  two  of  the  fiducials  (on  opposite  sides 
of  the  plate)  and  the  perpendicular  dropped  to  it  from  the  third  fiducial. 
Master  fiducial  coordinates  were  then  derived  by  computing  distances  from 
the  principal  point  to  each  of  the  measured  fiducial  images.  The  values 
used  for  preprocessing  were  based  on  the  mean  values  computed  from  measure¬ 
ments  on  ten  frames  of  Strip  3.  Although  the  master  fiducial  coordinates 
derived  from  Strip  3  were  applicable  to  Strip  4,  it  was  necessary  to  ad¬ 
just  two  of  the  values  by  approximately  0.100  mm  each  for  successful  pre¬ 
processing  of  Strip  5.  This  may  be  an  indication  of  significant  film 
distortion  in  Strip  5,  relative  to  Strips  3  and  4. 

Initial  aerotriangulation  of  each  strip  was  made  assuming  a 
0.010  mm  standard  deviation  for  x  and  y  on  the  plate.  The  0.010  mm  value 
was  chosen  as  a  starting  estimate,  although  not  a  very  realistic  one  for 
the  quality  of  input  imposed  on  the  triangulation.  As  expected,  the  re¬ 
sults  of  all  three  strips  indicated  that  the  plate  coordinate  weights 
should  be  relaxed  by  a  factor  of  approximately  four.  The  strips  were  re¬ 
run  assuming  0.040  mm  standard  deviations  in  the  plate  coordinates,  and 
results  indicated  that  consistent  weighting  was  used. 

3.6  Limitations  of  This  Stereo  Radar  Analysis 

Before  concluding  this  report,  a  few  words  should  be  said  about 
the  limitations  of  this  stereo  radar  analysis. 

Firstly,  the  analysis  was  conducted  using  very  selective  mater¬ 
ial.  Of  all  the  radar  imagery  acquired  for  the  test,  only  a  small  por¬ 
tion,  the  best  of  it,  was  selected  for  analysis.  This,  of  course,  was 
done  purposely  to  find  out  what  results  could  be  obtained  with  reasonably 
decent  imagery,  and  also  with  .the  expectation  that  future  imagery  would 
be  as  good,  if  not  better,  than  this  selected  imagery. 

Secondly,  for  reasons  given  earlier,  the  selected  test  area  was 
very  small  and  lacking  in  mountainous  terrain.  The  smallness  of  the  area 
limited  the  number  of  test  points  available  for  selection.  Only  49  points 
were  selected  for  Stereo  Radar  Model  I,  and  30  points  for  Model  II, 
hardly  enough  from  which  to  draw  conclusive  results,  but  nevertheless 
sufficient  to  give  useful  measures  of  accuracy.  The  lack  of  mountainous 
terrain  in  the  test  area  precludes  drawing  any  conclusions  from  this 
analysis  about  stereo  radar  accuracy  in  this  kind  of  terrain. 

Thirdly,  there  was  very  little  collateral  data  supporting  the 
radar  imagery  and  aerial  photography,  thus  requiring  adoption  of  a  special 
approach  to  data  reduction.  As  a  consequence,  the  results  are  directly 
related  to  the  particular  approach  used,  and  should  not  be  extrapolated, 
without  caution,  to  other  possible  stereo  radar  reduction  techniques. 
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Had  collateral  data,  such  as  range  and  time  marks,  been  furnished,  there 
is  really  no  telling  what  improvements  in  accuracy,  if  any  at  all,  would 
have  occurred. 

Fourthly,  the  photogrammetrically  triangulated  test  point 
positions,  used  as  the  standard  of  comparison  for  measuring  the  radar 
accuracy,  were  not  without  error  themselves,  inhibiting  measurement  of 
the  true  accuracy  of  the  radar. 

Finally,  the  specific  test  points  selected  for  analysis  were 
well-defined  and  very  carefully  identified  in  the  imagery.  They  were 
certainly  not  representative  ol  the  myriad  of  image  points  that  would 
have  to  be  used  for  preparation  of  a  topographic  map. 

In  spite  of  these  limitations,  this  analysis  still  gives  good 
insight  into  the  accuracy  of  stereo  radar.  Quite  clearly,  it  shows  that 
the  AN/APQ-102 (XA-2)  radar  has  the  geometrical  fidelity  necessary  for  the 
construction  of  topographic  maps  at  1:250,000  scale  and,  even  possibly, 
1:50,000  scale. 


CONCLUSIONS 


From  this  particular  stereo  radar  analysis,  the  following  con¬ 
clusions  are  drawn: 

(1)  The  AN/AFQ-102(XA-2)  radar  does  a  surprisingly  good  job 
of  locating  well— identified  ground  points  accurately 
using  stereo  radar  techniques.  In  general,  the  opposite- 
side  case  (Model  I)  provides  greater  accuracy  than  the 
same— side  case  (Model  II);  and  X  (the  across-track 
coordinate)  is  determined  with  greater  accuracy  than  Y 
(the  along-track  coordinate)  or  Z  (the  elevation). 
Specifically,  average  RMS  values  of  7.7  meters,  12.1  meters 
and  13.2  meters  in  X,  Y  and  Z,  respectively,  were  obtained 
for  Model  I;  and  average  RMS  values  of  9.5  meters, 

20.0  meters  and  16.7  meters  were  obtained  for  Model  II. 

(2)  The  geometric  fidelity  of  the  AN/APQ-102 (XA-2)  radar 
is  clearly  demonstrated  by  the  excellent  fit  of  the 
three  selected  radar  records  to  ground  control  provided 
by  photogrammetric  aerotriangulation .  Specifically, 
the  RMS  values  in  X  are  16.1  meters,  13.5  meters  and 
10.1  meters  in  Strips  3,  4  and  5,  respectively;  and  the 
RMS  values  in  Y  are  12.4  meters,  20.0  meters  and 

29.4  meters. 

(3)  The  horizontal  accuracy  provided  by  the  AN/APQ-102 (XA-2) 
radar,  using  stereo  radar  techniques,  is  adequate  for  all 
topographic  mapping  at  scales  of  1:50,000  and  1:250,000, 
except  perhaps  Class  A  mapping  at  1:50,000  seal  a.  Even 
for  this  one  possible  exception,  Model  I  (opposite-side) 
stereo  radar  may  provide  sufficient  accuracy. 

(4)  Vertical  accuracy  is  really  the  deciding  factor  with 
respect  to  applicability  of  the  AN/APQ-102 (XA-2)  radar 
to  topographic  mapping.  For  mapping  at  1:250,000  scale, 
the  Model  I  (opposite-side  case)  vertical  accuracy  seems 
sufficient  for  Class  A  mapping  with  50-meter  contours, 
or  for  Clasu  B  mapping  with  25-metcr  contours.  For 
1:50,000  scale,  the  Model  I  vertical  accuracy  seems 
sufficient,  only  under  the  best  test  conditions,  for 
Class  B  mapping  with  20-meter  contours,  or  for  Class  C-l 
mapping  with  10-meter  contours. 

The  conclusions  just  stated.,  are  based  upon  actual  results  of 
the  specific  tests  made  in  this  investigation.  It  is  conceivable  that 
with  better  imagery,  with  more  collateral  data,  and  with  more  accurate 
positions  for  the  standard  of  comparison,  the  results  could  be  even  better. 
On  the  other  hand,  for  a  larger  test  area,  with  mountainous  terrain,  the 
results  could  be  worse.  Therefore,  it  remains  to  be  seen  what  accuracy 
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can  be  obtained  over  an  extended  area,  with  varied  relief,  using  high- 
quality  imagery  with  all  possible  collateral  data. 

one  8eneral  conclusion  might  be  drawn  from  this  study,  it 
t^e  (XA-2)  radar,  using  stereo  radar  techniques,  par- 

t  cu  ar  y  with  opposite-side  configuration,  yields  enough  accuracy  to 
make  it  a  very  attractive  contender  for  the  job  of  all-weather,  day-night 
topographic  mapping  at  scales  of  1:250,000  and  even  1:50,000.  Of  course, 
more  tests  under  better  controlled  conditions  will  have  to  be  made  in  order 
to  arrxve  at  any  really  conclusive  results.  In  the  meantime,  the  results 
of  this  test,  qualified  as  they  are,  clearly  point  to  the  potential  of 
this  stereo  radar  technique  for  topographic  mapping. 
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5. 


RP COMMENDATIONS 


The  following  two  recommendations  are  made: 

(1)  More  extensive  tests  with  real  data  should  be  made, 
under  the  conditions  originally  anticipated  for  this 
investigation,  and  using  better  ground  control,  in 
order  to  obtain  more  conclusive  results. 

(2)  Theoretical  analyses,  using  simulated  radar  data, 
should  be  made  for  the  purpose  of  testing  the  effects 
of  various  errors  in  the  stereo  radar  model.  Such 
analyses  would  provide  insight  into  the  stereo  radar 
reduction  technique  with  respect  to  what  is  and  what 
isn't  important. 
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APPENDIX 


MATHEMATICAL  FORMULATION 


DETERMINATION  OF  RADAR  AIR  STATION  POSITION  AND  SLANT  RANGE 
FOR  A  GIVEN  CONTROL  POINT 

Let  C]  and  C2  be  two  camera  stations  initiating  and  termi¬ 
nating  a  line  segment  which  defines  the  flight  path,  Figure  A-l .  The 
radar  air  station,  R,  is  located  at  the  foot  of  the  zero-Doppler  line, 
which  is  normal  to  the  flight  path  and  passes  through  the  control 
point  I. 


jd  C  2  (X  ,  Y  ,  Z  ) 

/  L.  2  '“*2  ^  2 


r(xr’yr’zr2>> 

Flight  Path  - 


ci(VVzcr 


Z* 


> 


X 

-*• 


zero-Doppler  Line 


I  (control  point) 

(Xj.Yj.Zj) 


Figure  A-l 


The  direction  numbers  of  the  flight  path  are: 


ot  =  X  -  X_ 

U2  Uj 


3  =  Yc  -  Y 


y  zc  zr 

U2  Uj 


A-l 


Since  the  zero-Doppler  line  is  normal  to  the  flight  path,  it 
follows  that: 


a(xR-xx)  +  0(yr-y:)  +  yCZjj-Zj)  =  o. 


Also,  since  the  flight  path  is  a  straight  line. 


<XR  -  V  'f  (Yi<  -  Yc  > 


and 


<ZR  -  ZC,>  '  I  (YR  "  V  ' 

These  three  equations  are  solved  for  XD ,  YD  and  Z_.,  the  coordinates  of 

K  K  K 

the  radar  air  station.  The  solution  is: 


aS 


\  „ 2  4.  Q2  J.  +  X, 


Y„  = 


a  +  3  + 

3S 


R  a2  +  62  +  y2  Yci 


Zn  = 


YS 


R  a2  +  32  +  y2  + 


in  which 


S  =  a(X  -  X  )  +  3(Y  -  Y  )  +  Y(Zt  -  Z  ) 

Once  the  position  of  the  radar  air  station  is  known,  the 
slant  range,  SR^,  to  the  control  point  I  can  be  computed- 


SRj  =  [ (XR  -  xl)^  +  (Yr  -  Yj) 2  +  (ZR  -  Zj)2 ] 


2^2 


A-2 


2. 


SLANT  RANGE  CALIBRAT ION 


R:  =  r:  cos  6  +  t  sin  G. 

For  a  least  squares  solution,  the  condition  equation  for 
control  point  I  is: 

I  ASR  +  BSR^rIcos  (i  +  tjSin  0)  +  CsR(rjCos  0  +  tjSin  0)2-SRi  = 

In  this  equation,  SRj  Is  the  "observed"  quantity;  and  A  ,  B  ,  c 
and  0  are  the  unknown  parameters.  ^ 

This  condition  equation  Is  linearized  by  Taylor's  expansion, 
and  placed  in  the  matrix  form: 

AI  VI  +  A  +  e°  =  0 


in  which 


and  e 


V 


l  is  the  residual  in  SR^., 


B.  = 


3F,  3F, 


3asr  3bsr  3csf 


A  = 


6A 


SR 


6B, 


;r 


6C 


SR 


60 


3FI  9Fl 

90 


,  the  vector  of  corrections  to  estimates 
of  the  unknown  quantities. 


is  Fj  evaluated  using  approximations  for  the  unknown  parameters, 


The  partial  derivatives  are  evaluated  as  follows: 


3  SR, 


=  -1 


3F, 


3A 


=  1 


SR 


3Fi 

-  tj  cos  0  +  tj  sin  0 
jR 


3ri  , 

—  =  (rx  cos  9  +  t  sin  6) ' 
SR 


3Fi 

30 


=  BSR(tI  cos  9  ~  r  sin  0)  + 


+  2CSR^ri  cos  9  +  fci  sin  0)  (tj  cos  0  -  r 


sin  0) 


For  n  control  points,  the  set  of  linearized  condition  equations 
for  the  least  squares  solution  is: 

AV+BA+£°=0 

in  which 

A  =  I 

n 


vr 

~Br 

V2 

,  B  = 

b2 

,  and  = 

C2 

• 

• 

• 

• 

• 

• 

• 

• 

• 

V 

B 

n 

_  n 

L_  n  _ 

The  normal  equations  are: 

btwb  a  +  BTWe°  =  0  , 

in  which 

W  =  [aoa1]-1  =  I 

n  ’ 

since  A  =  I  ,  and  0  =  In;  i.e.,  the  computed  slant  ranges  are  assumed 
independent  and  equal  in  precision. 

The  normal  equations  are  solved  for  A: 

A  =  -  (BTWB)_1  BTtvE° 
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The  solution  is  iterated  with  the  parameter  approximations 
being  updated  by  the  A  corrections  for  each  iteration.  When  A  ■  0  or 
less  than  a  predetermined  convergence  limit,  the  solution  is  completed. 

The  residuals,  V,  are  then  computed  to  determine  how  well  the 
polynomial  fits  the  given  data. 

Coefficients  can  be  evaluated  for  each  radar  record,  or  the 

radar  records  can  be  combined  to  determine  common  values  for  B  and 

SR 

CSR.  In  this  case,  the  solution  incorporates  the  data  from  all  radar 
records,  and  an  expanded  set  of  condition  equations  is  used. 
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3. 


TIME  CORRELATION 


The  'time"  polynomials  for  point  I  are: 


\  ’  Ax  +  Bx  Ti  +  cx  V 


\  ■  AV  +  BV  h  *  CY  V 


\  *  AZ  +  BZ  TI  +  CZ  V 


in  which  XR  ,  Yr  and  ZR  are  the  radar  air  station  coordinates  for  the 
Ith  ground  point;  and  the  plate  coordinate,  Tj. ,  is  computed  as  follows: 


TI  ”  tl  cos  ®  ~  rp  sin  9  • 

The  coefficients  Ax>  Bx,  Cx>  Ay,  By,  Cy,  Az>  Bz>  and  Cz  are 
evaluated  using  the  method  of  least  squares.  For  each  control  point,  1, 
the  three  condition  equations  are: 


\  ■  Ax  +  Bx  TI  +  cx  V  -  Y  -  0 
\  ■  S  +  N  TI  +  CY  TT 


-  A„  + 


E  T 
Z  I 


C  T  ‘ 
Z  I 


-  Z. 


=  0 


These  condition  equations  are  already  in  linear  form  with 
respec :  to  the  unknown  coefficients.  Furthermore,  with  T  fixed,  and 
Xr^,  Yr^  and  ZR^  as  the  "observed"  quantities,  the  condition  equations 

are  independent,  and  thus  can  be  treated  separately. 


The  normal  equations  are  then: 


n  ETj  ET  2 

1 — ■  — 

*x 

EX 

RI 

ETj  ETj2  ETj 3 

Bx 

= 

ETI  \ 

ETx2  ETi3  ET  4 

cx 

^  _ 

These  normal  equations  are  solved  for  A^  Bx  and  Cx>  Similar  solutions 
can  be  obtained  for  the  Y-equations  and  the  Z-equations. 
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2)  Equation  of  plane  perpendicular 


to  Jj  and  through  I. 


where  oil,  3 1  ,  Yi  are  direction  numbers  for  Ji  at  the 
time  of  Ri. 


Two  similar  equations  can  be  written  for  J2 ,  resulting  in  a  total  of  4 
equations  ip  3  unknowns  (X^ ,  ,  z  ). 

Fi,  *  SRi,  -  |(xrxR1)2  +  <W2  +  (zi'V2ih ' 0 

°i i  '  “1<xrxR1)  +  <¥I-YRl)  +  y'<zi-zRi)  -  o 

fi2  ■  SRI2  -  [<xrV2  +  <VV2  +  (W2)15  -  0 

°i2  =  “2(xrV  +  62  (W  +  Y2(zrzR2)  -  o  . 


In  the  4  derived  equations,  a„d  SR^  are  functions  of  image  coordi- 

<3  t"  Q  C>  D  nn  J  H  .  1.1  *  • 


nates  ^  and  R^.  ^  ;  that  is, 


SRl!  =  ASRl  +  BSRl  \  +  CSRl  Rl' 
SRI,  =  ASR2  +  BSR2  RI2  +  CSR2  «!*.• 


Ri  *  *  ZRi  *  ^R2  *  YR2  an<^  ^r2  are  functions  of  image  coordinates 

T  and  T  ;  that  is, 

1  i  2 


Xr1  ■  **.  +  BX,  h,  +  CX, 

V  ■  **,  +  \  TI,  +  s, 


+  B 


Zi 


+  C 


Zi 


A-1I 


X 

Y 

Z 


r2 

r2 

r2 


V  +  BX2  ti2  +  Cx2 
=  \2  +  by2  ti2  +  cy2 


Image  coordinates  R  ,  T  R  and  Tt  are  considered  as 

1 1  J- 1  i  2  I2 

observations.  The  unknown  quantities  are  X^  ,  Yj  and  Z  . 

The  linearized  equations  in  the  general  matrix  form  are; 


Ai  h +  bi  Ai +  -  °  • 


where 

Aj  is  a  4x4  coefficient  matrix  of  partial  derivatives  of  the 
condition  equations  with  respect  to  the  observations. 

vj  is  a  4*1  column  vector  of  corrections  to  the  unadjusted 
observations. 

Bj  is  a  4X3  matrix  of  partial  derivatives  of  the  condition 
equations  with  respect  to  the  unknown  parameters. 

Aj  is  a  3*1  column  vector  of  corrections  to  the  parameter 
approximations. 

is  a  4*1  column  vector  of  constants  which  are  the  condition 
equations  evaluated  for  the  approximate  parameters. 

Now, 
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are : 


fAXT  " 

~  — 

F° 

I 

Ii 

G° 

ayt 

and  eJ = 

Ii 

F° 

I2 

azt 

G° 

L  i 

I2 

—  — > 

For  the  AT  matrix,  the  partial  derivatives  of  FT  and  FT 
1  i  1  -t  2 


3Ft  3Ft  3SR 

1 1  _  I  i  •  _ li 


t  -vT 

3SRt 

3Rt 

Ii 

1 1 

Ii 

3Ft 

3FT 

3SRt 

I2 

I2 

I2 

3Rt 

3SRt 

‘  3R 

I2 

I2 

I2 

3Ft 

3FT 

3X„  3Ft 

Ii 

_  Ii  . 

Rl  +  .  11 

3Ti, 

‘  3xr, 

3TI,  mr, 

3  F  _ 

3Ft 

3Xn  3Ft 

I2 

_  I2  . 

R2  .  I2 

3T 

J-2 

3V 

3ti2  8yr2 

3Y„ 

3Ft 

3Zn 

Rl  + 

Ii 

Ri 

3T 

3Z„ 

3Tt 

» 

Ii 

Ri 

Ii 

3Yn 

3Ft 

3Z„ 

r2  + 

I2 

.  Rz 

3Tt 

3Z_, 

3Tt 

> 

I2 

r2 

I2 

from  which 


3F 

3RT  '  2<ASRl  +  BSR,  RI,  +  CSR,  Rl!>(BSRl  +  2CSR,  RI C > 


3F, 

wf  -  2<asr2  +  bsr2  ri2  +  csr2  ri2>(bsr2  +  2csr2  ri2> 
I  2 


3F, 


X  - 


3Tt  SRt 
il  II 


<JF. 


■JL  - 


3T  SRt 

I  2  I  2 


(Xi-xEi)(bXi+2cXi  tIi)+(Vyr1)(Bv1+2cy1  V  + 


+  (Zi-ZRi)(Bz+2CZi  TXi) 


(xrV(Bx2+2cx2  V+<W<by2+2cy2  V  + 

'  <ZrZR2)(BZ2+2CZ2  V 


A-H 


Similarly,  the  partial  derivatices  of  G  and  G  are: 

X1  I2 

3G  3G 

_ Ll  —  Li.  _  o 

3R,  3R  U  * 

j.  i  I2 

and 


3°I 

IT^  -  -  ( V2CX,  v  -ei<Bv,+2cv,  V  -Yi  <BZl+2cZl  v 


3G. 


I2 


3t^  '  -  “><V2cx2  V  -6=  <by2+2cy2  V  nr*  <Bz2+2cz2  ti2> 


The  required  partial  derivatives  for  the  Bj  matrix  are: 


(xrxRi) 

9fi2 

<W 

3Xj. 

srt 

il 

3XI  ' 

srt 

12 

9FI, 

<W 

3fi2 

.  <VV 

3vi 

SRI, 

’  3Y 

srt 
i  2 

3Fi, 

.  (zrV 

3fi2 

(zi-zr2> 

3Z 

srt 

J-l 

'  3ZI 

SRi2 

3Gi, 

3xi 

=  Oil 

3gi2 
’  9x 

=  a2 

3Gi, 

3vi 

=  3i 

3gi2 

’  3Y 

=  g2 

3Gi, 

3Zj 

=  Yi 

3GI2 
'  3zi 

=  Y2 
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The  least  squares  solution  is  as  follows: 


The  weight  matrix  W  is 


wi  -  <Ai  °i  Ah" 


where 


a 


I 


o 


in  which  the  a2,s  are  the  estimated  variances  of  the  image  coordinates. 


The  correction  vector  Aj  is  computed  from  the  matrix  equation 
AXt 


AI  = 


AY, 


AZ, 


-  -<BI  wi  v~'  «i  q 


The  solution  is  iterated  with  the  parameter  approximations  being 
updated  by  the  A  corrections  for  each  iteration.  When  A  =  0  or  is  less 
than  a  predetermined  convergence  limit,  the  solution  is  completed. 

A  statistical  analysis  of  the  solution  can  be  performed  by 
making  the  following  computations; 


The  vector  Vj  is  computed  from 
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VI  =  -(°1  Ar  V1)~l  (Bj.  Aj  +  ej) 


VI  ai  VI  is  comPuted  directly,  or  from  =  e°T  y  e° 

The  a  posteriori  variance  of  unit  weight  02Oj  can  then  be  found  from 


T  -i 

Vt  °t  v  T 

J2  III 


’I 


M  -  P 


where  M  is  the  number  of  condition  equations  and  P  is 
the  number  of  unknown  parameters. 


The  covariance  matrix  for  the  unknown  parameters  can  be  estimated  from 
cov  Aj  =  (b]:  Wx  Bj)'1  . 
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